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ABSTRACT 

We study the evolution of the intracluster medium (ICM) with a uniformly analyzed sample of 70 
galaxy clusters spanning 0.18 < z < 1.24 and observed with Chandra. We find that X-ray luminosity 
and ICM mass at a fixed temperature evolve with redshift in a manner inconsistent with cither the 
standard self-similar model of cluster formation or a model that assumes no evolution of cluster 
structure. Both luminosity and ICM mass evolve more slowly than the self-similar prediction, i.e., 
clusters have lower luminosity and ICM mass at fixed emission-weighted temperature than expected at 
higher redshifts. We find that evolution in these two observables can be modeled by a simple evolution 
in the cluster gas mass fraction, evolving as (1 + z )~°- 39±013 when measured using core-subtracted 
observables. Excluding cluster cores from measurements results in evolution more consistent with 
the self-similar model than when the entire cluster is used, indicating that the fraction of clusters 
with cool cores increases with time, or that cool cores become more developed over time in those 
clusters that have them; this is supported by direct study of the redshift dependence of central surface 
brightness, which increases in scatter and magnitude at low redshift. We find that isophotal size- 
temperature relations evolve differently according to which isophote is used, indicating that the central 
and outer regions of cluster ICM evolve differently. We show that constraints on the evolution of the 
gas fraction and isophotal size-temperature relations constraints can be combined to measure cluster 
distances, and thus to constrain cosmological parameters in a way complementary to other techniques. 
Scatter in scaling relations is considerably reduced by using cither core-subtracted quantities or three- 
parameter relations including the central surface brightness; in addition, there are indications that 
scatter decreases at higher redshift, suggesting that merging is not the dominant source of cluster 
structural variation. Our results provide constraints for simulations attempting to model cluster 
physics, indicate some difficulties for cosmological studies that assume constant cluster gas fractions, 
and point toward other potentially more robust uses of clusters for cosmological applications. 
Subject headings: galaxies: clusters: general — X-rays: galaxies: clusters — intergalactic medium 



1. INTRODUCTION 

Scaling relations among bulk properties of galaxy clus- 
ters provide a powerful means to test models of the large- 
scale structure and evolution of the universe. These 
correlations among properties such as X-ray luminosity, 
intracluster medium (ICM) mass, mean ICM temper- 
ature, and cluster virial mass reflect gravitational and 
non-gravitational processes involved in the formation of 
structure in an expanding universe. Scaling relations also 
provide the means to readily estimate masses of clusters 
from much more easily measured properties such as lu- 
minosity, an essential component of cosmological stud- 
ies that use X-ray observations to determine the redshift 
evolution of the cluster mass function. 

Simple models of cluster formation via gravitational 
collapse predict particular forms for the redsh ift evolu- 
tion of cluster scaling relations (Kaiser 19861. Adding 
additional cluster physics such as radiative cooling of 
the ICM, and energy injection by active galactic nuclei 
(AGN), superno vae, and star formation, modifies these 
predictions (e.g.,|Cavaliere et al 
Muanwong et al 



studies of scaling relation evolution are required to prop- 
erly constrain models of cluster evolution and to under- 
stand the effects of non-gravitational processes on the 
scaling relations that will be used to study cosmology. 
X-ray studies of the ICM are complementary to studies 
of the evolution of the cluster galaxy population (e.g. 



de Propris et al. 
strain the overal 



1999| |Lin et al.||2006[ ), helping to con- 
evolution of cluster baryons and their 



distribution in various forms within clusters. 

Several studies of X-ray scaling relatio n evolution have 
been carried out in recent years (e.g., 
[20021 |Ettori et al. 
Maughan et al.|2006 
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120071) , but no clear consensus has emerged. In this paper 
we will address scaling relation evolution using a system- 
atic analysis of a Chandra sample of 70 clusters covering 
0.18 < z < 1.24, the largest sample yet used for this 
purpose. 

Our study addresses two difficulties which may affect 
scaling evolution measurements. The first arises from 
the fact that radiative cooling of the ICM leads to the 
development of cool, dense (and hence very luminous) 
cores in many clusters; these relatively small cores bias 
cluster measurements such as X-ray temperature and lu- 
minosity to an extent that they are not representative of 
the overall cluster structure. This introduces significant 
scatter into scaling relations; indeed, there is evidence 
that cool core clusters, which are traditionally regarded 
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as "relaxed," actually exhibit greater structural variation 
than non-cool core clusters, which are oft en thought to 



have recently undergone major mergers (O'Hara et al 
2006|) . Studies of scaling relations commonly attempt to 



"correct" for the impact of cool cores on cluster prop- 
erties by one of several methods, such as simply leaving 
clusters wit h evidence for strong co ol cores out of the 
sample (e.g., Arnaud fc Evrard|1999 1, or excising central 
region s within a fixed metric radius (e.g., Morandi et al.| 
2007| or a fraction of the virial radius (e.g., |Mau ghan 
et al. 12007 1, and perhaps "correcting" measured luminos- 



ity by some factor determined from a m odel of the clus 



ter su rface brightness distribution (e.g., Vikhlinin et al 



2002 1 . In this paper we measure temperatures with and 



without cores defined as fractions of the virial radius, and 
we also measure luminosities with and without the same 
core. By using relations both with and without core sub- 
tracted quantities, we can examine the effects that core 
development has on cluster scaling relation slopes and 
evolution. 

The other issue usually faced by scaling relation stud- 
ies is the use of scaling relation slopes and normaliza- 
tions from low-redshift studies carried out with differ- 
ent instruments. The relatively small fields of view of 
Chandra and XMM-Newton make measurements of local 
samples quite challenging with those instruments; hence, 
studies using older X-ray instruments are used as refer- 
ences for z = relations. Unfortunately, differences in 
spectral and imaging results among X-ray instruments 
are well established, making such approaches subject to 
instrument-related systematics; indeed, even the same 
instrument has produced results differing by the author, 
as calibrations change and varying reduction and anal- 
ysis methods are adopted. By using a large sample (70 
clusters) , we can avoid the use of outside references for 
scaling relation parameters or the direct inclusion of data 
from other samples, in favor of a single, homogeneously 
analyzed sample. While this approach is not entirely 
new — for example, Branchesi et al. (20071 studied evo- 
lution using their own 17 cluster sample both with and 
without the inclusion of data from other studies; and 
|Morandi et al.| (|2007[ ) studied a homogeneously reduced 
24 cluster sample — the size of our sample leads to signif- 
icantly smaller uncertainties on scaling relation parame- 
ters than have otherwise been obtained. 

In § [2] we provide a brief overview of scaling relations 
and their predicted evolution, and in §[3] we explain our 
data reduction and measurement procedures. We test 
for scaling relation evolution with respect to expecta- 
tions from the self-similar theory and from a scenario of 
no evolution in cluster parameters in §§[4] and [5] respec- 
tively, and provide an explanation for observed evolution 
in scaling relations via a simple evolution in the gas mass 
fraction §|6] In §[7] we examine the evolution of isophotal 
size, and discuss the implications for studying cosmology 
using size measurements, and in §[8] we discuss the effec- 
tiveness of two different methods of reducing the scatter 
in measured scaling relations. In § [9] we compare our 
results to previous observations and simulation results, 
and discuss some implications of our findings. Finally, 
we list our conclusions in § IP] 
We adopt th e WMA P + LRG ACDM cosmology from 



Spergel et al. (20071, which combines the third year 
WMAP data with results from the SDSS luminous red 



galaxy survey ( Eisenstein et al.|[2005 ) to give H = 70.9 
km s^ 1 Mpc" 1 , il M = 0.266, and ii A = 0.734. All uncer- 
tainties are 68% confidence, unless specified otherwise. 

2. SCALING RELATION BACKGROUND 



The self-similar model (e.g., Kaiser|1986 1 describes for- 
mation of clusters via gravitational collapse of overdense 
regions in an expanding universe. In this model the ICM 
is heated by this gravitational collapse and the resulting 
shock heating, but no so-called non-gravitational heating 
is assumed. As a result, clusters scale self-similarly, i.e., 
they scale only because of changes in their physical size 
at fixed mass due to density variation as the universe 
expands. With the assumptions of spherical symmetry, 
hydrostatic equilibrium, a constant gas fraction, and X- 
ray emission dominated by thermal bremsstrahlung, this 
leads to X-ray luminosity Lx and ICM mass M g scaling 
with ICM temperature Tx and redshift as 



L x cx T x E(z), 



M g tx T x 



3/2 



E(zY 



(1) 



(2) 



where E(z) is the ratio of the Hubble parameter at red- 
shift z to its present value. In a flat cosmology with 
matter density f2 m , E(z) has the form: 



E{z) = H(z)/H = [n m (i + z) 3 + 1 - n m ] 



1/2 



(3) 



Predicting scaling laws for the isophotal size (i.e., the 
physical size of the region corresponding t o th e angular 
size of a particular X-ray isophote; see § 3.5 1 requires 
additional assumptions about the ICM mass distribution. 
If the ICM distribution scaled self-similarly with mass, 
then isophotal size scales as 



i?j cx Tx /3 , 



(4) 



with no redshift evolution if the ICM density falls off as 
outside the core (i.e., /? = 2 



Mohr et al. 2000). 



Observational studies have found that scaling relations 
for all three of these observables (Lx, M g , and Ri) in 
fact have a stronger dependence on t emperature than 
predicted by self-similar models (e.g., |Edge fc Stewart 



T99Tj |Markevitch|r9 98 ; Mo hr fc Evrard[lQ97[|Mohr et al 
1999[ ) . Explanations for this and other evidence of non- 
gravitational processes, such as the presence o f entropy 
ramps in t he central regions of clusters (e.g., Ponman 
et al.|2003 1 , typically involve additional energy injection 
by active galact ic nuclei (AGN), supernovae, and star 



formation (e.g ., Bialek et al.| 
I Borgani et aTplX)^ |McCarth 



2001 Bower et al. 2001 



try et al.||2004| |Kay et al| 
|2007[ ); radiative cooling of the ICM, which leads to the 
formation of cool, dens e cores in ma ny clusters; and non- 
radiative cooling (e.g., Bryan 2000). 

It is important to note that there are multiple ways 
to define radii for measuring cluster parameters, which 
result in different predicted redshift evolution for scal- 
ing relations. The expressions given above are correct 
for observables (Lx and M g ) measured within regions 
corresponding to fixed overdensities relative to the crit- 
ical density. This is appropriate for our strategy in this 
paper, in which we choose to measure cluster proper- 
ties within virial regions defined by local relations, and 
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then test for consistency with the evolution scenarios de- 
scribed below. Another commonly used form for the 



rcdshift evolution of scaling relations (e.g., 


Ettori et al. 


2004b 


Branchesi et al. 


2007 


Morandi et a 


. 2007p uses 



in a spherical collapse model. These densities have their 
own redshift evolution, leading to additional factors in 
the scaling relation evolution equations. In cither case, 
it is common to parametrize additional rcdshift evolution 
beyond the self-similar predictions in terms of a simple 
power law with redshift, i.e., proportional to (l+z) raised 
to some power. 

In this paper we discuss two models for cluster evolu- 
tion. The first is "self-similar evolution," in which cluster 
observables scale as would be expected given purely grav- 
itational influence as discussed above, i.e., Lx oc E(z) 
and M g oc E(z) . The other is what we will refer to 
as "no evolution," meaning that cluster parameters, in- 
cluding virial radii, do not scale at all as the universe 
expands. 



> 
v 




Fig. 1. — Measured emission-weighted mean temperature Tx 
plotted versus redshift for the clusters in our sample. 



3. DATA REDUCTION 
3.1. The Cluster Sample 

The data are drawn from the Chandra archive. The 
lower redshift limit of z ~ 0.2 reflects the difficulty in 
measuring cluster parameters out to at least r 2 5oo for 
clusters closer than this, given the small Chandra field of 
view. The cluster sample is listed in Table [6] with the ID 
number of the Chandra observation used for each cluster. 

Having been largely developed through cluster selec- 
tion in archival Einstein IPC and ROSAT PSPC ob- 
servations, our sample is essentially X-ray flux limited. 
However, as the sample is not derived from a single ho- 
mogeneous survey at a fixed flux threshold, it might be 
expected to include systematically more luminous (i.e., 
more massive) systems at high redshifts. In Figure [I] we 
plot the emission-weighted mean tem peratures for our 
sample (measured as described in £3.3 below) versus red- 



shift. Our sample spans a consistent range of Txover the 
full redshift range. 



3.2. X-ray Data Reduction 

The data reduction is carried out using the standard 
Chandra analysis software CIAO, version 3.3, with CALDB 
version 3.2.1, and the spectral fitting package xspec, ver- 
sion 11.3.1. We generate new level 2 events files from the 
level 1 files obtained from the Chandra archive, so that 
all observations are reduced in a uniform manner. The 
following reduction procedure is applied to each cluster. 

Light curves are extracted for back-illuminated chips 
5 and 7 individually, and for front-illuminated chips 0-3 
and 6 combined. Light curves are extracted and binned 
in time using the recommended criteria for each chip. 4 
Flares are excluded using the CIAO task "lc_CLEAN" 
based on the median value of the light curve. The expo- 
sure times after filtering are given in Table [6j 

Cosmic ray events arc identified with the CIAO tool 
"acis_RUN_HOTPIX" . A new level 1 events file is then 
generated using the latest gain file, and charge trans- 
fer inefficiency (CTI) and time-dependent gain variation 
corrections are applied as appropriate. Standard bad 
columns and hot pixels are excluded. Events with ASCA 
grades of 0, 2, 3, 4, and 6 are used. A level 2 events file is 
then created from the filtered level 1 events file. Where 
the observation was made in very faint (VF) mode, we 
carry out the extra background event flagging that this 
enables. 

We attempt to use background data from the actual 
data sets, extracting the background from regions well 
away from target cluster or other emission. For some 
clusters, however, emission fills most of the detector, and 
in these cases we extract the background spectrum from 
the Markevitch blank-sky data. 5 To account for small 
differences in the particle background between these sta- 
tistical backgrounds and each individual observation, the 
blank-sky sets' exposure times are scaled by the ra- 
tio of counts in the 7-12 keV energy band in the data 
and blank-sky observations. Before using either back- 
ground method point sour ces are identified by the iter- 
ative method described in Sanderson et al. ( |2005[ ) and 
checked by visual inspection, and then excluded. Even 
when emission-free regions are available, if the spectral fit 
is worse with the local background than with the blank- 
sky background, we use the latter. In total, we use the 
blank-sky backgrounds for 41 of the 70 clusters in our 
sample. 

3.3. Spectral Fitting 

Cluster spectra are extracted in regions with maxi- 
mum radius chosen by eye to be where the cluster emis- 
sion merges into the background; the center coordinates 
and radii of our extraction regions are given in Table [6] 
Choosing apertures based on the X-ray surface bright- 
ness distribution might result in smaller apertures rela- 
tive to the physical size of clusters that are cooler or lie at 
higher redshifts, and thus will tend to have observations 
with fewer total counts. However, Figure [2j which plots 
the ratio of the spectral extraction radius to r5oo for each 
cluster versus cluster mean temperature (left) and red- 
shift (right), suggests that this is not the case. The mean 
ratio of aperture radius to rsoo is 0.84±0.20 (RMS), with 
no apparent temperature or redshift dependence. 

4 http://cxc.harvard.edu/ciao/ 

5 http:/ /cxc. harvard.edu/contrib/maxim/acisbg/ 
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We generate weighted response matrix files (RMFs) 
using the CIAO tool mkacisrmf when the data allow; 
otherwise we use the older tool mkrmf. 

We fit to the cluster spectra a single-temperature 
APEC model with a component for galactic absorption. 
We use fit -/Vh values when they are reasonable (i.e., 
within a few standard deviations of the galactic value), 
and not p egged to zero; otherwise, w e fix Nn to the galac- 
tic value (Dic key fc Lockman||1990 1. In total, we fit JVh 
for 18 of the 70 clusters. We generally extract spectra in 
energy bands of 0.7-9 keV for ACIS-I, and 0.5-9 keV for 
ACIS-S. In a few cases we use an upper limit of 7 keV 
when there is clearly spurious, non-background emission 
above this value; in no case does this change the mea- 
sured temperature at g reater than the 1-2% level. We 
use the Cash statistic (|Cash 19791, which is preferable 



to x 2 when the number of photons is low. In our sample 
the use of the Cash statistic generally results in a best- 
fit temperature that is a few percent higher than that 
measured with \ 2 . 

We measure the core subtracted temperature Txcs by 
extracting spectra with the same maximum radius as de- 
scribed above, but excluding the inner 0.2r 50 o; the core 
subtracted temperature and the 0.2r5oo exclusion radius 
are measured iteratively until convergence. (Our defini- 
tion of rsoo is given in § |3.5|) For two clusters, ZwCl 



1356+6245 and CLG J0647+7015, the iteration does not 
converge to a reasonable value when the core is excluded, 
and so we do not measure core subtracted quantities for 
those two clusters. 

Our measured values for the temperature of the entire 
cluster, and for Txcs measured assuming self-similar evo- 
lution and assuming no evolution, are given in Table [6] 

3.4. Comparison with Published Temperatures 

Though calibrations continue to improve, measure- 
ments of the same cluster by different instruments, and 
by different methods with the same instruments, lead to 
temperature measurements that differ. To check the ac- 
curacy of our own temperature measurements, we com- 
pare our values to those obtained in two other recent 
Chandra studies. 



Balestra et al. (20071 (hereafter Ba07) studied 56 clus- 
ters over a redshift and temperature range similar to our 
own; our samples have 38 clusters in common. Our data 
reduction and spectral fitting processes differ from theirs 
in several small ways: Ba07 use local backgrounds ex- 
clusively, while we, as described above, use blank-sky 
backgrounds when local backgrounds are not possible or 
give worse spectral fits; they always fix the value of Nn 
to galactic, while we allow it to float when the value 
obtained thereby is reasonable; they use a spectral ex- 
traction band of 0.6-8 keV, versus our 0.5- or 0.7-9 keV; 
and they include a spectral component to compensate 
for Ir-M edge residuals, a correction that has been taken 
into account in the more recent calibration files which we 
use. Because clusters are not isothermal, the emission- 
weighted mean temperature is affected by the choice of 
energy band. Most importantly, Ba07 use spectral ex- 
traction regions determined via a method intended to 
maximize the S/N, which results in the use of extrac- 
tion radii up to a factor of two smaller than ours. Their 
resulting extraction regions have a clear redshift trend, 
with radii as small as ^0.3 rsoo at high redshift. 



Maughan et al. ( 2007 1 (hereafter Ma07) measured tem- 
peratures for 115 clusters, of which 53 are in common 
with our sample. Differences between our analyses in- 
clude their use of a 0.6-9.0 keV spectral fitting band; 
their fixing iVjj to the galactic values; and their use of 
blank-sky backgrounds in some cases where we use local 
backgrounds, plus an additional soft X-ray background 
component. Ma07 also use a different method for de- 
termining the spectral extraction region, measuring all 
spectra out to a radius of r§oo as determined from an it- 
erative procedure using a mass-lx relation, w here Yx is 



the produc t of the temperature and gas mass (Kravtsov 
eTil1[2006| 



To examine the difference between our temperatures 
and those of these two studies, we compare the error- 
weighted ratio of our temperatures to theirs. Overall, our 
temperatures are lower than those of Ba07 by a weighted 
average of (3 ± 1)%, and higher than those of Ma07 by 
(6 ± 1)%. To examine whether we can reproduce their 
values, we remeasured the temperatures of five clusters 
using methods similar to those of Ba07 and Ma07; i.e., 
we used their reported aperture radii, spectral extraction 
bands, and spectral models. We fixed Nr in all cases 
for this comparison, but did not change our choice of 
background strategies. As can be seen in Table [l] these 
changes resulted in generally higher temperatures when 
using the methods closer to those of Ba07, and generally 
lower temperatures when using methods closer to those 
of Ma07, thus at least partially explaining the sources 
of systematic differences between our measurements and 
those of these two papers. Note that this does not mean 
that our temperatures necessarily came to agree more 
closely with theirs; for MS 0451.6-0305, for example, our 
initial temperature was higher than that of Ba07, and 
these changes resulted in an even higher temperature. 

The overall hotter temperatures that we measure rel- 
ative to Ma07 may be attributable to variations in ICM 
temperature with radius. As shown in |3.3| our spectral 
extraction radii average (0.84 ± 0.20)r5oo, while Ma07 
uses uniform radii of t^qq. The ICM temperat ure gener- 
ally d ecreases with radius at these radii (e.g., |Vikhlinin| 



et al. 2005 Pratt et al.||2007| ), and so we would expect 
our measured temperatures to be systematically slightly 
higher than those of Ma07. 

However, the differences between our temperature 
measurements and those of the other two studies are not 
uniform; there arc dependences on temperature and, for 
Ba07, redshift. The left panel of Figure [3] shows the ra- 
tio of our temperatures to the literature values versus 
our temperature. In the case of the Ma07 comparison, 
the ratio is clearly greater at higher temperatures; a one- 
dimensional least-squares fit of a straight line shows that 
the ratio increases as (0.14^ 19) log Tx for Ba07, and as 
(0.22i°;°g)logT x for Ma07. The latter trend may again 
result from Ma07's choice of rsoo as an extraction ra- 
dius; extraction regions of hotter clusters may include 
more background-dominated area, leading to tempera- 
ture systematics as parts of the spectrum are dewcighted 
by background noise. 

The right panel of Figure [3] shows the dependence of 
temperature ratio on redshift. There is no evidence for a 
redshift dependence when comparing to Ma07; the ratio 
varies as (0.02 ± 0.09)2;. For Ba07, however, the ratio 
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Fig. 2. — The ratio of our spectral extraction radius to rsoo (defined as described in j ]3.5[ l for each cluster, plotted versus the measured 
non-core subtracted temperature (left) and versus redshift (right). Markers correspond to z < 0.4 (circles), 0.4 < z < 0.6 (squares), and 
z > 0.6 (triangles). 



TABLE 1 

Comparison of Temperature Measurements 









Balostra ct al. (2007) 


Maughan et al. (2007) 


Cluster 


z 


Our T x (kcV) 


They (kcV) Wc" (koV) 


They (kcV) Wc" (keV) 



MS 0451.6-0305 


0.54 


9.8 ± 0.8 


8 9 + 0-4 


10.5 ± 0.7 


6 7 +0 ' 6 

°' ' -0.5 


8.1 ± 0.4 


C1G J 1149+2223 


0.54 


9.8 ± 0.8 


12.9+^ 


9 9 +1 


8 4+ ' 9 

0,<l — 0.7 


a 7 + 0.9 


C1G .11120+2326 


0.56 


4 2+ - 6 


5.2 ± 0.5 


4 4+0-4 


, a + 0.4 
3.»_ .3 


3.2 ± 0.3 


C1G J1113-2615 


0.73 


, 7+0.6 
a - 1 -0.5 


- 7 +0.9 


5 +0 ' 9 


o q+0.9 
3-8-0.7 


3.1 ± 0.4 


RX J1317+2911 


0.81 


3 8+ 17 


4 5+ 1 - 4 


4 4+0-9 


2 0+ ' 7 


o q + 1.5 

°°~o.s 



a Our measurement of the cluster temperature using the same aperture and similar methods as the 
literature sources; sec text. 



varies as (— 0.31^Q gg)z, showing a clear negative depen- 
dence on redshift. This is almost certainly a result of 
Ba07's use of extraction regions that feature a trend to- 
ward smaller fractions of the virial radius at higher red- 
shift. 

The differences between our measured temperatures 
and those from the literature underscore the difficulties 
inherent in comparing cluster parameters measured us- 
ing differing instruments, instrumental calibrations, and 
methods. This calls into question the reliability of re- 
sults obtained from directly combining data from mul- 



tiple studies (e.g., |Branchesi et al. 2007J), and suggests 
that caution should be taken when comparing more pro- 
cessed results, such as the low-redshift slopes and nor- 
malizations often combined with new measurements of 



tion (e.g., Vikhlinin et al. 2002 


Ettori et al. 2004b 


Kotov 


& Vikhlinin||2005||. 



3.5. Imaging Analysis 

We extract X-ray images and use the spectral fit to ob- 
tain the conversion factor from counts to physical units 
in the rest frame 0.5-2 keV band. Because the flattening 
of statistical backgrounds using the exposure map gener- 
ated for a particular observation results in a spatially in- 
homogeneous background image, we fit a flat background 
to the regions outside of the cluster emission using the 
same technique used to determine the surface brightness 



profile, described below. The results of this fitting are 
checked by examining radial brightness profiles and via 
simple comparison of total counts in regions well outside 
of cluster emission. 

As our observations do not in general contain enough 
photons to do a deprojection analysis, particularly at 
high redshift, we fit the standa rd spherical (3 model ( Cav- 
aliere fc Fusco-Femiano|[l978 1 to the cluster emission: 



I(r) = h 



-i -3/3+1/2 



(5) 



with central brightness 1$, core radius R c , and power-law 
index (3. In what are traditionally considered "cool core" 
clusters, i.e., where there is a central emission excess due 
to the formation of a coo l dens e core, we fit a dou ble 



model (llkebe et al.||1996l 119991 iMohr et al.||1999|) with 



both components having the same center coordinates and 
index (3, so that the total surface brightness is the sum 
of the two, i.e., 



i=i 



R r 



-3/3+1/2 



(G) 



We fit these surface brightness profiles to the two- 
dimensional surface brightness images, and find the best 
fit and one a confidence intervals for each parameter us- 



O'Hara et al. 




Fig. 3. — The ratio of our measured cluster temperatures to publi shed temperatur es (vertical axis ) , plotted versus our t emperature (left) 
and versus redshift (right). Published temperatures are from |Maughan etTa l. (2007) (triangles) and Bal estra et al.| l |2007| ) (circles). 



ing the Cash statistic. In a few cases cases (A521, A1682, 
and A2744) there are prominent clumps or subclusters 
separate from the main body of the cluster, which are 
masked out before fitting. In two cases w e fix the value 
for (3 : A521, a multiply- merging cluster (Ferrari et al. 
20031, for which we find a somewhat stable value of 
(3 = 0.75, which we adopt over the values of (3 > 3 which 
are found by a full grid ding analysis; an d C1G J 1056 



0337, a merging system (Jee et al 



20051 for which we 



r 500 = 1.129 



r 2 500 = 0.501 



5 keV J 
5 keV J 



0.497 



E(z)- 1 Mpc, (7) 



0.503 



E(z)- 1 Mpc. 



(8) 



Note that by using virial radii obtained in this man- 
ner, we are implicitly testing the evolution of these local 
mass-temperature relations along with our other observ- 
ables. That is, our "self-similar evolution" scenario in- 
cludes evolution of the r$-Tx relations as written above, 
and the "no evolution" scenario includes no evolution 
(i.e., no E(z) factor) in the r^-Tx relations. 



find only very high values of f3, leading us to adopt the Ferrari et al 
canonical (3 — 0.67. In both cases we then measure 2 a 
uncertainties in the other fit parameters. The (3 model 
parameter fit results are listed in Table [7] The second, 
bright central component is, where used for a given clus- 
ter, listed as the second brightness and core radius com- 
ponents I2 and R Cy 2- 

We measure several different cluster observables, each 
of which — X-ray luminosity, ICM mass, isophotal size, 
and mean ICM temperature — derives from the under- 
lying cluster structure in a different way; by studying 
the evolution of multiple observables, we are examining 
the evolution of the ICM in multiple ways. Luminosity 
and ICM mass arc measured within two different virial 
radii r$ , which permits us to examine evolution on differ- 
ent scales within a cluster. We determine r§Qo and r25oo 
from the clu ster temperature using M$-Tx relations de- 
termined by Arnaud et al. (20051 using XMM-Newton 
observations of local galaxy clusters. We use their rela- 
tions for clusters with Tx > 3.5 keV: 



We measure the projected X-ray luminosity Lx in the 
rest frame 0.5-2 keV band from the images described 
above, within radii of r5oo and ^2500; we also measure 
core subtracted luminosities Lxcs by excising the pro- 
jected luminosity from the central 0.2r5oo- Luminos- 
ity measurements are centered on the cluster brightness 
peak, with the exception of A521, where we use the peak 
brightness of the main cluster, not the brighter infalling 
subcluster to the n orth of the cluster center (see, e.g., 
and C1G J1056-0337, where we use 



2006) 



the western brightness peak, which has been identi fied 
as the "cen tral" mass peak via weak lensing (e.g., Jee 
et al. 2005 ) . Given the small field of view of Chandra, 
the virial radii r$ often extend beyond the image bound- 



ary; furthermore, some observations are not deep enough 
that there is signal measurable out to a given r$. We 
thus establish for each cluster a maximum radius from 
the brightness peak at which either the detector edge is 
reached or the S/N falls close to unity; in a few cases 
the maximum radius is determined by the presence of 
other structure, as in the cases of ACO 2246 and C1G 
J1701+6414, which lie a small angular distance from one 
another in the same observation. Then, if the radius 
rg exceeds this established maximum radius for a given 
luminosity measurement, we do not carry out that mea- 
surement on that particular cluster; this is reflected in 
Tables [8] and [9] where luminosity measurements are not 
given in many cases. We include in the luminosity un- 
certainties contributions from the temperature used in 
calculating rg, as well as a uniform 10% background un- 
certainty. 

The X-ray luminosity within a given radius can be 
modeled analytically by an integral of the ICM density 
profile and X-ray emissivity out to that radius. We can 
therefore use a measurement of the actual luminosity to- 
gether with the measured (3 model parameters and the 
cluster temperature to find the central ICM density, and 
hence ICM mass via an integral of the density function 
to a g iven radius of interest; for details see |Mohr et al. 
(19991. We estimate uncertainties on M g by including 
the statistical uncertainties on the (3 model fit; a uniform 
10% background uncertainty in the luminosity measure- 
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ment; and temperature uncertainties in r$. The ICM 
mass measurement is not subject to the same maximum 
radius restriction as luminosity, as the luminosity within 
any given radius can be used to measure the central den- 
sity; while larger luminosity measurement radii are of 
course preferable, it is not necessary to measure the flux 
out to a given rs for an ICM mass measurement within 
that radius. Note that we do not similarly use the (3 
model to extrapolate luminosity measurements out to a 
radius of interest; this is because we prefer to directly 
use projected luminosities without assumptions as to the 
structure of the cluster, but ICM mass cannot similarly 
be measured without such assumptions. 

We measure the isophotal size Rj of a cluster by mea- 
suring the area Aj enclosed by an isophote /, and find- 
ing the effective radius given by Aj = nRj. For these 
measurements we use images that have been adaptively 
smoothed using the CIAO task csmooth. We include the 
10% background uncertainty in the Ri uncertainties by 
remeasuring at isophotes increased and decreased by the 
background uncertainty. In the 0.5-2 keV band we are 
using here, the conversion from X-ray counts to physical 
units varies slowly with cluster temperature, so we do not 
include temperature uncertainties in the isophotal size. 
We measure Ri at three different isophotes, 1.5 x 10 -13 , 
6 x 10 -14 , and 3 x 10~ 14 erg s _1 cm -2 arcmin" 2 (in the 
rest frame 0.5-2 keV imaging band), which, like using 
both rsoo and r25oo for the luminosity and ICM mass, lets 
us study evolution of Rj on different scales within a clus- 
ter. Clusters can "fall off" an isophotal size-temperature 
scaling relation when the isophote used approaches the 
peak surface brightness of the cluster; we therefore ex- 
clude clusters when their measured isophotal size is less 
than 0.2r 50 o, our adopted core exclusion region. 

3.6. Fitting Procedures 

For a given relation involving an observable O, we fit 
the form 



O = A 



6 keV 



(9) 



or, in log space, 



log O = log A + a log 



6 keV 



+ 7log(l + z). (10) 



That is, we fit a power-law temperature dependence a, 
power-law redshift dependence 7, and A, the normaliza- 
tion at zero redshift and temperature 6 keV. 

In this paper we use unweighted orthogonal fits, mean- 
ing that we minimize the sum of the square of the point- 
line orthogonal distances, i.e., the sum 



E 



logO, - [logA+alog(T x ,i/6) + 7 log(l + z t )} 



(1 + a 2 ) 



1/2 



(11) 

Note that the form for redshift evolution assumed here 
is evolution of the normalization only, and so there is no 
factor of 7 in the denominator. We determine 1 a uncer- 
tainties via bootstrap sampling; the best-fit value given 
in this paper is the mode of a histogram constructed from 



the bootstrapping results, and the 1 a confidence inter- 
val is constructed in the usual manner so as to contain 
68.3% of the counts around this mode. We also give here 
the RMS scatter in the vertical dimension (e.g., in Lx 
in the Lx~Tx relation) for the best-fit parameters; this 
one-dimensional scatter is a more intuitively understand- 
able quantity than the orthogonal scatter, as it reflects 
the scatter in an observable (Lx, M g , Ri) at a given 
temperature. We refer to this as the intrinsic scatter 
((7i n t), as the measurement uncertainties are generally 
much smaller than the total scatter in these relations 
(e.g., |Q'Hara et al.|[2006| . 

The question of which fitting method is "best" is still 
open, and rests to a large extent on whether one property 
(such as Tx) is considered more fundamental than the 
other (such as Lx)', this often seems implicit in discus- 
sions of Lx~Tx, Mg-Tx, and other relations, and would 
imply that a one-dimensional least-squares fit , with tem- 
perature (the lowest scatter mass estimator) as the in- 
dependent variable, might be appropriate. But if both 
observables are considered to be linked via another prop- 
erty of the system (such as cluster mass) , then a orthogo- 
nal minimization fit, which treats both variables equally, 
may be more appropriate; we take this view, and so adopt 
orthogonal fitting in this paper. 

Fits of mock scaling relations using both the orthogo- 
nal fitting method and an ordinary least-squares (OLS) 
fit support this decision. A difficulty that arises in such 
tests is that assumptions must be made regarding the 
scatter in mock relations; e.g., if only scatter in the y di- 
rection is generated, then an OLS fit will doubtle ss give 
better result s than an orthogonal fit. For example, |Lopes| 
et al. ( 2006} make the claim that orthogonal regression 



produces more accurate measurements of scaling relation 
slopes than the bisector method (discussed below) , based 
on their own tests using mock data sets; however, as these 
data sets were generated using orthogonal scatter, such 
a result is entirely expected. Because of the difficulty 
in defining "correct" scatter, we test scenarios in which 
only scatter in the y direction is used, and in which equal 
scatter in both the x and y directions is used. That is, 
we generate a random x value, use an assumed scaling 
relation to find y, and then shift the values using normal 
random deviates in the y direction only, or in both the x 
and y directions. Note that using equal x and y scatter 
is not the same as using orthogonal scatter, and so an or- 
thogonal relation should not be a priori assumed to give 
the correct result in such a case. In our testing we do 
not assume measurement uncertainties, but fit an intrin- 
sic scatter in the y direction in the OLS fits so that the 
reduced x 2 value is equal to unity. Again, in real scaling 
relations the scatter is generally dominated by intrinsic 
scatter, so this is a reasonable approach. 

The results of our tests clearly indicate that the OLS 
method is a less robust approach than the orthogonal 
method. For example, when using only y direction scat- 
ter of 0.05 (i.e., the random deviates have a standard 
deviation in log 10 space of 0.05), the orthogonal method 
gives a result that is 2% (~ler) high while the OLS 
method gives the correct result; but when using equal 
scatter of 0.05 in x and y, the orthogonal method gives 
the correct slope, while the OLS method gives a result 
that is ~10% (~2cr) too low. The results get worse for 
OLS more rapidly than for orthogonal fitting; e.g., scat- 
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Fig. 4. — Projected X-ray luminosity within r2soo (left) and r5oo(right), with non-core subtracted (top) and core subtracted (bottom) 
quantities, plotted versus temperature. These quantities are measured assuming self-similar evolution. Luminosity values are scaled to 
z = using the best-fit redshift scaling for each relation, and the best-fit slope is plotted for each relation. Markers correspond to z < 0.4 
(circles), 0.4 < z < 0.6 (squares), and z > 0.6 (triangles). 
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Fig. 5. — Same as Figure [4] but for M 9 -Tx relations. 
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ter of 0.15 in y only gives an orthogonal slope that is 
16% (~2.5cr) high, but scatter of 0.10 in both dimen- 
sions gives an OLS result that is 51% (~8.5c) too low. 
The results are very similar when true orthogonal scat- 
ter is used, rather than random, but on average equal, 
scatter in each dimension. 

Again, the exact origin of scaling relation scatter is 
unknown, so it is difficult to declare a "correct" way 
of testing fitting methods. There is undoubtedly some 
measurement scatter, however, and so scaling relations 
certainly have at least some scatter in both dimensions. 
For this reason, as well as the physical arguments given 
above, we adopt the orthogonal fit as our chosen method 
for this paper. 

Besides orthogonal fitting, another approach that 
treats the two variables equally is the bisector method, in 
which OLS fits are done with each of the two variables 
as independent and dependent (i.e., y as a function of 
x, and function of y) , and the final result bisects 

the two individual fits. This is not appropriate for our 
work, because we fit observables as a function of both 
temperature and redshift, and it is unclear how the bi- 
sector method can be extended into three dimensions. 
Orthogonal fitting is clearly defined in any number of 
dimensions; i.e., it seeks the shortest point-line distance 
in two dimensions, the shortest point-plane distance in 
three dimensions, and so forth. Also, each individual 
OLS fit in the bisector method is subject to the great 
dependence on the form of scatter as discussed above, 
and so the bisector method's utility for studying scaling 
relations is likewise questionable. 

4. TESTS OF THE SELF-SIMILAR EVOLUTION SCENARIO 

We now examine the evolution of scaling relations 
while assuming self-similar evolution, as discussed in the 
introduction. That is, we assume that r$ scales as E(z)~ 1 
when measuring Lx and M g , and when determining the 
core subtraction radius for Txcs and Zoccs- Our values 
for Lx and M g , measured using the non-core subtracted 
temperature, are given in Table [8] We then test whether 
scaling relations evolve in a manner consistent with self- 
similar evolution. 

4.1. Scaling Relations 

The Lx~Tx and M g -Tx scaling relations are plotted 
in Figures [4] and [5] respectively. In these figures the 
observables are scaled to z = using the best-fit scal- 
ing relations. One qualitative feature of note is that 
the scatter is clearly smaller in the Lxcs~Tx relations 
than in their non-core subtracted counterparts; a simi- 
lar, though smaller effect is visible in the ICM mass re- 
lations. This difference in scatter arises from biases in 
both temperature and the other observable in each rela- 
tion i nduced as a result of cool core-related phenomena 



relation and Lx,500~2x relation have slopes of 2.75_q' 2 q 

and 2.35±g:||, respectively, which are significantly higher 
than the self-similar expectation a = 2, as h as been gen- 
erally observed in low-redshift samples (e.g., |Markevitch 



Fabian et al.|1994| |Markevitch||1998| |Q'Hara eTal 



2006[ ). Another interesting feature is the shallowness of 
the -£xcs,2500~2xcs relation compared to the non-core 
subtracted Lx,2500~Tx relation. Best-fit scaling relation 
parameters are given in Table [2] 

Studies of scaling relation evolution commonly fix the 
slopes to values measured from local samples, and fit 
only for an evolution factor. Because we are fitting all 
parameters simultaneously, we need to compare our mea- 
sured slopes to those of local samples. Our Lx,25oo _ ?x 



craliy c 

19981; note that using luminosities from a fixed energy 
band as done here (rest frame 0.5-2 keV) gives a some- 
what lower slope than the more commonl y used bolomet - 
ric luminosities, as shown by, e.g., |Zhang et al" (20071. 
For the M gj 25oo _ Tx relation we find a = 1.82 ± .08, in 
good a greem ent with a — 1.91 ± 0.16 found by |EttorT] 



et al. (20021 using BeppoSAX data and a bisector fit; 
for the Mg t500 -Tx relation we find a = 1.74 ± 0. 09, in 
fair agreeme nt with a — 1.98 ± 0.11 measured by Mohr 
et al.|(|1999|) using ROSAT PSPC images and a mixture 



of Einstein, Ginga, and ASCA temperatures, with an un- 
weighted orthogonal fit. Both of these are significantly 
higher than the self-similar expectation a = 1.5. 

In all cases, the scaling relations with core subtracted 
quantities have shallower slopes than the standard rela- 
tions. Remarkably, the core subtracted relations have 
slopes consistent with the self-similar expectation to 
within lcr, the sole exception being M 9i 2500~2xcs, which 
is consistent to within 2a. 

4.2. Evolution with Redshift 

Figure [6] shows the ratio of observables (Lx and M g ) 
to the self-similar expectation, plotted versus redshift. 
That is, the vertical axis is the ratio of the observed 
value to the self-similar prediction using the appropriate 
fit in Tableland the cluster temperature and redshift. 
The horizontal line in each plot therefore marks the self- 
similar expectation. Plotting in this way shows devia- 
tions from the self-similar redshift evolution prediction 
as a redshift dependence of the ratio Oj/Ofit = 0; we 
also plot the best-fit value of 7 for each relation, show- 
ing how the normalization of each scaling relation in fact 
evolves. 

For each scaling relation, Table [2] includes the percent 
significance by which each relation differs from zero, i.e., 
the significance of its deviation from the self-similar pre- 
diction. Because the distributions of 7 are not in gen- 
eral normal, this significance is determined using binned 
data to measure the probability density at 7 = 0, and 
integrating to the same probability density on the other 
side of the peak value. Because we use binned data to 
estimate this parameter, it can be determined most pre- 
cisely when 7 is significantly different from zero; hence, 
we quote only at 1% precision for values less than 99%. 

All luminosity- and ICM mass temperature scaling re- 
lations have 7 < at greater than the lcr level. There is 
clearly an overall tendency for relations to evolve more 
slowly than expected from the self-similar prediction, i.e., 
7 < 0. We can combine multiple probabilities by assum- 
ing independence of the scaling relations; though all of 
the measured properties are of course linked to some ex- 
tent by their dependence on the underlying ICM struc- 
ture, X-ray luminosity and ICM mass depend on that 
structure in very different ways, and the two virial radii 
which we use probe two rather different regions of the 
cluster (i.e., rsoo comes close to looking at the cluster as 
a whole, while r25oo measures a much smaller fraction 
that is more dependent on core structure and evolution) . 
Combining the results for all four core subtracted rela- 
tions by multiplying the given probabilities of consistency 
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TABLE 2 



Fit Parameters Assuming Self-Similar Evolution 


Relation 


a 


Core Subtracted Relations 
A a 7 


Diff. from 

(%) b 




-kxcs,250o-Txcs 

-^XCS ,500~ TxCS 
A^g,2500~^XCS 
^g^OO-^XCS 


1 00+ 023 

2 26+ ' 29 

z -^°-0.33 
1 fi Q+0.09 

- L - DC> -o.os 
1.56 ± 0.10 


2.00+. ;^E44 
3.02lJ;| 4 E44 
2.57lg;J?E13 
7.94i -° 6 5 E13 


-0.86±°;1° 

-o -«+ - 2 ° 

-0 24+ ' 20 

u '^*— 0.18 


95- 
69- 
90- 
74- 


0.28 ± 0.05 
0.21+"™ 

— 0.07 

00+ 05 
0.09 ± 0.04 


Relation 


Non-Core Subtracted Relations 
a A a 7 


Diff. from 
(%) b 


C^int C 


£x,2500-Tx 

Afg,2500 — 
M gi 50Q—Tx 


„ -,r+U.29 

z - ,o -0.26 
9 ,,-+0.33 
'• JO -0.24 

1.82 ± 0.08 
1.74 ± 0.09 


3.24+JJ;^E44 
6.03t 2 ' 22 E44 
2.69t°i|E13 
8.32^° ;^E13 


-1 50 +u - 42 
-1 90+ 117 

-0 ^ + 0-17 

-0 45+ ' 18 


99.4- 
90- 

99.4- 
98- 


60+ u u8 
39+ 12 
0.14 ± 0.02 
0.13 ± 0.02 



a In units of Lq for Lx-Tx relations, Mg for M g — Tx relations. b Significance level at 
which 7 differs from zero, as determined by bootstrap sampling and refitting; the sign 
indicates whether 7 is positive (+) or negative ( — ). c Intrinsic scatter in Lx or M g at fixed 
temperature, expressed in base e. 




Fig. 6. — Ratio of measured observable (luminosity or ICM mass) to the best-fit observable-temperature scaling relation, plotted versus 
redshift. These measurements assume self-similar evolution. The horizontal line (O/Ogt = 1) corresponds to no evolution beyond the 
assumed self-similar evolution, i.e., 7 = in our notation. The dashed and dotted lines correspond to the best fit and 1<t boundaries on 7 
for each relation. 



with zero gives a combined probability of < 0.1% that 
all four relations are consistent with the self-similar evo- 
lution scenario, ruling out pure self-similar evolution at 
greater than 3er confidence. The same relations with non- 
core subtracted quantities have an even smaller proba- 
bility (i.e., <C 0.1%) of consistency with zero. 

We draw your attention to the z > 0.8 clusters in 
our sample because of the special leverage they have on 
our evolution results. Examination of Figure [6] suggests 



no qualitative difference in the high-rcdshift population 
when compared to lower-redshift clusters. For these clus- 
ters to bias our results toward more negative evolution, 
it would require systematically selecting underluminous 
clusters, which is the opposite of what is expected. 

The relations involving core subtracted quantities have 
more positive evolution than those involving non-core 
subtracted quantities. This could indicate a decrease in 
clusters with cool cores at higher rcdshifts, which is ex- 
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pected in the scenario wherein clusters form cool cores 
over time in the absence of major merging events. The 
evolution of the cool core frac tion remains relatively un- 
explored; Bauer et al. (2005) found no evolution in the 
cool core fraction up to z ~ 0.4 using spatially resolved 
spectral analysis, but suc h an analysis is difficult to carry 
out at higher redshifts. Vikhlinin et al. (20061 used a 
measurement of the "cuspiness" of the surface brightness 
distribution to count cool cores in a sample of clusters 
at z > 0.5, and found a fourfold decrease in the cool 
core fraction from z=0 to z=0.5, which might support 
the concept of cool cores indicating a "relaxed cluster" 
that has not undergone recent major mergers. This con- 
cept is being increasingly challenged, however, by results 
from simulations that ascribe the presence or lack of a 
cool core to aspects of cluster form ation history such as 
pre heating (McCarthy et al.||2004 ) or early major merg- 



ers (|fjurns et al.| |2007p , and observational evidence that 
on-cool 



cool core and non-cool core cluster populations differ in 
characteristi cs beyond the ir morphological state ( |Q'Hara| 
et al. 20061. Burns et al. (20071 specifically studied the 
redshift evolution of the cool core fraction, and find no 
change in the fraction up to z ~ 1 in simulations that 
successfully reproduce other aspects of cluster and core 
structure. Our results here may support the classical no- 
tion of cool cores evo lving over time, in support of the 
Vikhlinin et al. (20061 results. Alternatively, a constant 
cool core fraction could still produce an apparent neg- 
ative evolution in scaling relation normalization simply 
because cool cores in those clusters that do have them 
will tend to grow over time; such a result was reported 
in simulations by Kay et al. (20071. We further discuss 
possible evolution in scatter in § [§] 



4.3. Summary of Self-Similar Evolution Results 

X-ray luminosity and ICM mass at fixed temperature 
evolve more slowly than expected from the self-similar 
evolution model. This conclusion is supported by sig- 
nificant (> la) negative evolutions in all Lx~ and M g - 
Tx scaling relations, and by combined constraints us- 
ing multiple core subtracted or non-core subtracted re- 
lations that rule out self-similar evolution at > 99.9% 
confidence. The less negative evolution of the core sub- 
tracted relations suggests that the cool core fraction de- 
creases with redshift, that cool cores grow over time, or 
a combination of the two. 



5. TESTS OF THE NO EVOLUTION SCENARIO 

We now examine the evolution of scaling relations 
while assuming no evolution, i.e., we assume no scaling 
in r$ when measuring Lx and M g , and when determin- 
ing the core subtraction radius for Txcs and Lxcs- Our 
values for Lx and M g measured using the non-core sub- 
tracted temperature are given in Table [9] We do not 
measure Lx.500 or -Lxcs, sou m this scenario, as only a 
handful of clusters have observations of sufficient expo- 
sure time and angular extent that we can measure out to 
the non-evolved r§oo- 
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Fig. 7. — Projected X-ray luminosity within r2soo (left) and 
r 500 (right) , with non-core subtracted (top) and core subtracted 
(bottom) quantities, plotted versus temperature. These quantities 
are measured assuming no evolution. Luminosity values are scaled 
to z = using the best-fit redshift scaling for each relation, and the 
best-fit slope is plotted for each relation. Markers correspond to 
z < 0.4 (circles), 0.4 < z < 0.6 (squares), and z > 0.6 (triangles). 

5.1. Scaling Relations and Their Evolution 

The Lx~Tx and Lx~M g relations are plotted in Fig- 
ures|7]and|8] respectively. As in the self-similar evolution 
case, the slope of the luminosity-temperature relation de- 
creases significantly when core-subtracted quantities are 
used, and the scatter likewise decreases for both the lumi- 
nosity and the ICM mass relations. We give the best-fit 
scaling relation parameters from this scenario in Table [3] 
and plot the redshift evolution of the scaling relations in 
Figure [9j 

The measured slopes and normalizations in this sce- 
nario are consistent with those measured in §[4] including 
the tendency for core subtracted relations to have shal- 
lower slopes than non-core subtracted relations. Also in 
common between the two scenarios is the tendency for 
core subtracted relations to have more positive evolution 
than non-core subtracted relations. 

Single non-core subtracted relations arc generally con- 
sistent with negative evolution, and core subtracted re- 
lations are generally consistent with positive evolution. 
Combining all three core subtracted relations gives a 
combined consistency with 7 = (i.e., with the predic- 
tions of the no evolution scenario) of 1%; for the non-core 
subtracted relations, the value is 8%. 

5.2. Summary of No Evolution Scenario Results 

The core subtracted scaling relations rule out the "no 
evolution" scenario at 99% confidence; non-core sub- 
tracted relations give less certain results. As in the self- 
similar evolution scenario, the core subtracted relations 
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TABLE 3 

Fit Parameters Assuming No Evolution 



Relation 


a 


Core Subtracted Relations 
A a 7 


Diff. from 

(%) b 


int C 


-kxcs,2500— Txcs 

^g,2500— ^XCS 


1 84+°' 18 
1.57 ± 0.07 
1.52 ± 0.08 


1.74l°;=|E44 
2.40t° ;^E13 
7.59j:S'5tE13 


n rc+0.37 
°- 56 -0.35 

20 +0 20 
10+ 019 


88+ 
89+ 
41 + 


24 +0 04 
0.00 ± 0.00 
0.08 ± 0.04 


R.clation 


Non-Core Subtracted Relations 
a A a 7 


Diff. from 

(%) b 


Cint C 


^x,2500— ^x 

^g,2500 — 


„ --+U.34 
z - ,o -0.25 

1.78 ± 0.08 
1.74 ± 0.10 


2-8lt^;™E44 
2.57+°;l 8 E13 
8.13+n"^E13 


-0.25 ± 0.56 
-0 10+ 017 
-0 25+ ' 19 


36- 
43- 
78- 


0.59 ± 0.09 
0.13 ± 0.02 
0.14 ± 0.02 



a In units of Lq for Lx—Tx relations, Mq for M g -Tx relations. 13 Significance level at which 
7 differs from zero, as determined by bootstrap sampling and refitting; the sign indicates 
whether 7 is positive (+) or negative ( — ). c Intrinsic scatter in Lx or M g at fixed temperature, 
expressed in base e. 
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Fig. 8. — Same as Figure [7] but for M 9 -Tx relations. 

have slopes that are consistent with self-similar expecta- 
tions, and evolution that is more rapid than the corre- 
sponding non-core subtracted relations. Together with 
the results from the self- similar evolution tests, these 
findings indicate that cluster scaling relations do evolve, 
but they evolve less rapidly than the self-similar expec- 
tation. 

6. TESTING EVOLUTION OF THE ICM FRACTION 

One simple model for the evolution of cluster parame- 
ters such as Lx and M g is a simple evolution of the gas 
mass fraction f g , i.e., the ratio of the ICM mass to the 
total mass (baryons + dark matter) of a cluster. It is 
sometimes assumed in cosmological studies using clus- 
ters that /„ is constant with redshift if clusters are se- 



lected appropriately (e.g., |Rines et al.||1999 [ [Allen et al. 



2004 1 , but this assumption is difficult to test because of 



degeneracies between f g measurements and cosmological 
parameters. Simulations disagree on the baryon fraction 
evolution, with some clai ming to see a negative evolu- 
tion ( e.g.,|Kay et al.|2007[ ), while others find no evolution 



(e.g., |Cram et al ||2007p . 

We can test whether our data are consistent with an 
evolution in f g by directly combining measured values of 
7 for individual scaling relations. X-ray luminosity varies 
proportional to the square of the ICM density, and ICM 
mass is directly proportional to the ICM density. Be- 
cause we are working in log space, this means that we 
combine jM g with 7l x /2. We use the core subtracted 
relations for this test because these relations are presum- 
ably less biased by cluster structural changes in the core, 
and therefore more sensitive to more global changes in 
the gas fraction. 
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Fig. 9. — Ratio of measured observable (luminosity or ICM mass) to the best-fit observable-temperature scaling relation, plotted versus 
redshift. These measurements assume no evolution. The horizontal line (O/O^t = 1) corresponds to no evolution, i.e., 7 = in our 
notation. The dashed and dotted lines correspond to the best fit and 1<t boundaries on 7 for each relation. 
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Fig. 10. — Constraints on the evolution of / 9 for the self-similar evolution (left) and no evolution (right) scenarios. Open triangles 
are from the fit to the Lxcs 2500"^XCS relation (with the values halved, as discussed in the text), open squares are from M g 2500-Txcs 
relation, and filled circles are the normalized product of the two. The best-fit to the combined relations gives y* = —0.39 ± 0.13 in the 

self-similar evolution scenario, and yt = 0.25^q 'Jj in the no evolution scenario. 



First we examine the ixcs,2500- and Af ffj2 500-Txcs 
relations measured in the self-similar evolution sce- 



nario. The left panel of Figure 10 shows histograms for 
the values of 7 resulting from the bootstrap fitting of 
the £xcs,2500 - an d Mg t 25oo-Tx relations (triangles and 
squares, respectively; the values of 7 for luminosity have 
been divided by 2 as explained above); the vertical axis 
has been scaled so that the values represent the proba- 
bility of 7 falling in each bin. The circles are the product 
of the two individual distributions, renormalized so that 
the total probability is unity. The data give a best fit 
value of — —0.39 ± 0.13; the data are inconsistent 
with jf = (i.e., a constant gas fraction) at the 99.1% 
level. 

The right panel of Figure [10] shows data calculated 
in the same way, but in the no evolution scenario. In 



this scenario we find the best-fit combined scaling to be 
7/ = 0.25^0^, and inconsistent with yj g — at the 
98% level. 

Our results are consistent with the evolution in Lx 
and M g originating from a simple evolution in gas mass 
fraction. While such consistency does not prove this sce- 
nario, it is encouraging to note that the values of jm and 
7l x /2 are quite similar in both scenarios, and evolution 
in f g thus provides a consistent explanation for the evolu- 
tion of these two different physical quantities. The most 
probable value 7/ ~ —0.4 in the self-similar evolution 
scenario suggests a decrease of ~25% in f g between red- 
shifts and 1. which would bias distance measurements 



that assume constant f a at the ~17% level (gU cx /, 



2 /3. 



e.g., |Rines et al ||1999[ ). 

Note that we have measured the evolution of f g specif- 
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Fig. 11. — Isophotal size— temperature relations for non-core subtracted (top) and core subtracted (bottom) temperature; the isophote 
used decreases from left to right. Size values are scaled to z = using the best-fit redshift scaling for each relation, and the best-fit slope 
is plotted for each relation. Markers vary by redshift as in Figure [4] 
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Fig. 12. — Ratio of measured isophotal size to the best-fit size— temperature scaling relation, plotted versus redshift. These measurements 
assume no evolution. The horizontal line {Ri / Ri fit = 1) corresponds to no evolution, i.e., 7 = in our notation. The dashed and dotted 
lines correspond to the best fit and lcr boundaries on 7 for each relation. 



ically within the radius r^soo- We do not attempt a sim- 
ilar measurement at rsoo because of a lack of luminosity 
measurements at that radius in the no evolution scenario, 
and the very large uncertainties on the £xcs,500 - Ixcs re- 
lation in the self-similar evolution scenario. There is in 
both scenarios and in both core subtracted and non-core 
subtracted relations a tendency for M 9 .5oo to evolve more 
slowly than M„ 2500 (though only at the 0.5-ler level); this 
is consistent with observations and simulations which 
find that the evolution in f g decreases with increasing 
radius , with evolution nearing zero at t he virial radius 
(e.g., |Sadat et al.|2005|rEttori et al.|2006[ ). This suggests 
that the distance biases associated with the assumption 
of constant gas fraction would be less severe if the X-ray 
data were deep enough to allow the measurements to be 



made at or beyond the virial radius, which may be pos- 
sible which future observatories such as Constellation- X. 

7. EVOLUTION OF ISOPHOTAL SIZE 

We now examine the evolution of isophotal size- 
temperature scaling relations. This is done separately 
from the previous "self-similar evolution" and "no evo- 
lution" because as discussed in § [2j for clusters that are 
described by a j3 model with (3 = |, the two scenar- 
ios give the same result (Mohr et al. 20001. That is, 
if the ICM is distributed as r _z beyond the core, the 
size-temperature relation does not evolve with redshift. 
While this potentially makes the size-temperature rela- 
tion useful as a means to study the evolution of the an- 
gular diameter distance, and hence as a tool for studying 
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TABLE 4 

Fit Parameters For Isophotal Size Relations 



Core Subtracted Relations 



Diff. from 



Relation 


a 


A (Mpc) 


7 


(%)" 




fl 1.5xlO- 13 ~ Txcs 

fl 6x 10 -i4-7xcs 


70+ 07 

u - ' u — 0.06 
0.08 


0.31 ± 0.02 
0.46 ± 0.03 


-0.05lg-.il 


42+ 
33- 


0.13 ± 0.02 
0.16 ± 0.02 


R 3 x r o - 1 4 -Txc s 


0.65 ± 0.10 


n fi o+0.05 
uo ' : >_0.04 


-0.26 ± 0.18 


87- 


n , 9 + 0.03 
U - 1Z -0.02 


Relation 


Non-Core Subtracted Relations 
a A (Mpc) 7 


Diff. from 

(%r 


_ b 




0.81 ± 0.07 


0.32 ± 0.01 


-0-03±«;i» 


13- 


0.13 ± 0.01 




0.81 ± 0.09 


0.47 ± 0.03 


-0.16 ± 0.16 


68- 


0.16 ± 0.02 




n 74+0. os 
u. (-r_ n n7 


60+ 06 


-0.22 ± 0.17 


79- 


0.14 ± 0.02 



a Significance level at which 7 differs from zero, as determined by bootstrap sampling and 
refitting; the sign indicates whether 7 is positive (+) or negative ( — ). b Intrinsic scatter in 
Ri at fixed temperature, expressed in base e. 



cosmology, it makes it less useful for constraining the 
evolution of the ICM and cluster structure as we have 
done with luminosity and ICM mass relations. 

7.1. Scaling Relations and Their Evolution 

Size temperature scaling relations are shown in Fig- 
ure [Tl] as with the previous scaling relation plots, these 
havenad the measured redshift evolution projected out. 
Best-fit scaling relation parameters are given in Table [4] 
The slopes of the relations using core subtracted tem- 
peratures are consistent with the theoretical value a = | 
( Mohr et al.||2000 l, and the relations with non-core sub- 
tracted temperature are somewhat higher. Our fit slope 
for the R 3xW -i4-T x relation is 0.74±°;^, which differs 
significantly from the value a = 0.93 ± 0.11 found by 
Mohr et al'.| fl2000[ ) using ROSAT PSPC images and lit- 
erature values for Tx . 

Redshift evolution of the isophotal size relations is 
shown in Figure [12] For the fits to the entire sample, 
the isophotal size relations show little or no evolution in 
the isophote closest to the core, and a trend toward more 
negative evolution as the isophote used decreases, i.e., as 
one examines the cluster at distances further from the 



|6 that the evolution in ix an d M„ 



core. 

Having shown in 
with respect to the self-similar expectation can be mod- 
eled by a simple evolution in the gas fraction, we can 
check for consistency of that evolution with the isopho- 
tal size results. The brightness at a given cluster radius 
r is related to the gas fraction f g as I(r) cx fg, and so it 
can be shown that for a cluster described by a spherical 
j3 model the meas ured isophotal size scales with I(r) as 
Ri x /( r ) i/(6/3-i) ( [Mohr et al.|2000| . Thus we expect 



R^f'^K 



(12) 



which, for the standard value of ft = 4 (e.g., Jones & 



size should scale as fg^ 3 . For our self-similar evolution 
measurement of jf = —0.39 ± 0.13, this would pre- 
dict Rj x (1 + z)~°' 26±0 - 09 , in good agreement with 
the directly measured evolution of 7 = —0.26 ± 0.18 



Forman 1984 Mohr et al. 19991, means that isophotal 



in the R? 



-Txcs relation, and of 7 



-0.05 



+0.15 
-0.18 



in the i?6xio^ 14_ 2xcs relation. More positive evolution 
at higher isophotes may be an indication of structural 



changes as clusters evolve and the density profiles of clus- 
ters become more peaked at lower redshift. 
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Fig. 13. — Measured angular diameter distance dj± plotted 
versus redshift for the ^3x10-14 sample. The solid line shows 
(Ia{z) for the best-fit measured cosmology Qm = 0.02, Qa = 0.31, 
and the dotted line shows dj^(z) for our adopted cosmology Hj^f = 
0.266, Q A = 0.734. 



7.2. Prospects for Cosmology Using Isophotal Size 

As mentioned above, the predicted non-evolution of Ri 
with redshift makes these size measurements a promising 
source of angular diameter distances, which can be used 
to constrain cosmological parameters. Such an under- 
taking is beyond the scope of this paper, but we sketch 
here the basic ideas underlying such a measurement. 

If isophotal size indeed evolves in a manner predictable 
by the evolution in f g , then one can use a measured an- 
gular isophotal size 9j, together with a physical isopho- 
tal size Ri for the same cluster predicted from a scal- 
ing relation, to determine the angular diameter distance, 
d,A = Ri /Oi ; this can then be used to measure the cosmo- 
logical parameters which determine oIa(z). As a test, we 
use 9 1 measured from our i?3 X io- 14 sample, and use the 
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best-fit slope and normalization found for the i? 3xl0 -i4- 
Txcs relation to predict Ri(Tx,z). Because we have 
found evolution in Lx and M g which suggests evolution 
in f g , we adopt the best-fit f g evolution 7/ s = —0.39 and 
its consequent isophotal size evolution 7^ = —0.26 in 
the size-temperature relation, as discussed above. Un- 
certainties in (La are a combination of the temperature 
uncertainty and the measured intrinsic scatter in the 
^3xio- 14 " T xcs relation. 

Note that this is not simply an independent cosmo- 
logical test. This is a consistency test where (1) an in- 
put cosmological model is assumed, (2) f g evolution is 
measured using the evolution-sensitive Lx and M g scal- 
ing relations, and (3) that evolution is adopted in using 
isophotal sizes to derive an output cosmology. Because 
the cosmological dependencies of each scaling relation 
differ, the input and output cosmologies will only agree 
for the correct model. 

Figure [13] shows the angular diameter distance versus 
redshift, with our input cosmology and the output best- 
fit cosmology. As can be seen, these data do not reach 
rcdshifts high enough to place tight constraints on cos- 
mology. Figure [14] shows confidence intervals for the den- 
sity parameters sTj^ and J7a (we fix Hq to our assumed 
value of 70.9 km s _1 Mpc -1 ). The uncertainties on both 
parameters are quite large; fully marginalized constraints 
are Q M = 0.02±g;g|, fl A = 0.31±g;f|. We do successfully 
recover our input cosmology within the 1 a region. 

This combination of the use of Lx~Tx and M g -Tx rela- 
tions to constrain the evolution of the ICM, and Rj-Tx 
relations to measure distances is an approach that de- 
serves further attention. As X-ray surveys that include 
spectroscopic temperature measurements push to higher 
rcdshifts, the use of isophotal sizes to measure angular 
diameter distances as demonstrated here should provide 
a new source of cosmological measurements, complemen- 
tary to other cluster methods and to CMB and supernova 
constraints. 
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Fig. 14. — Constraints on Qm and Qa from fitting the angu- 
lar diameter distances determined from isophotal sizes. The thick 
and thin contours mark the boundaries of the 1 and la confidence 
regions, respectively. The circle denotes the best fit fij^f = 0.02, 
£7a = 0.31, and the cross marks our adopted cosmology for this 
paper, Cl M = 0.266, fi A = 0.734. 



8. SCATTER IN SCALING RELATIONS 

This paper has focused on the evolution of the nor- 
malization of observable-temperature scaling relations. 
Here we briefly discuss the scatter about those scaling 
relations, i.e., the variation in the ICM distribution from 
cluster to cluster at fixed temperature. Understanding 
the precise origins of scatter helps both in understanding 
cluster physics such as cool core development and merger 
effects, and in understanding sources of uncertainty in 
cosmological studies that use observables such as X-ray 
luminosity an d temperature as prox ies for cluster mass. 
As shown by O'Hara et al. (20061, the cluster central 
surface brightness Iq is strongly correlated with central 
cooling time and reflects the core structure of clusters. In 
this section we examine the use of Iq to reduce scatter in 
scaling relations, and to examine the redshift evolution 
of cluster structure. 

8.1. Reducing Scatter: Two Approaches 

As shown in previous sections, the total scatter in scal- 
ing relations generally decreases when core-subtracted 
quantities are used, reflecting the separation in cool core 
and non-cool core pop ulations that is observed in most 
sc aling relations 



■g 



Fabian et al. 1994 Markevitch 



1998| [McCarth y et al.|2004"||O'Mara et al.|2006p . JO'Hara 



[et*aT] ( P006p demonstrated that central surface brightness 
Iq can be used as a proxy for cool core "strength" in a 
three parameter (0-Tx~Iq) scaling relation, reducing the 
scatter in scaling relations that is introduced by biases to 
both the temperature and to the other observable (Lx, 
M g , Rj) in the relation. With the data presented here we 
can compare the three-parameter approach to the use of 
core subtracted quantities, to determine whether either 
method results in lower scatter than the other. 

Rather than usin g the f3 model values for Iq, as in 
O'Hara et al. (2006 1, we estimate Iq by simply averaging 
the surface brightness within 0.05r5oo of the brightness 
peak. Since our intention is to use Iq to parametrize 
the development of cool cores, this method is likely to 
give more accurate results than the surface brightness 
fitting which, even when a double [3 model is used, may 
not accurately reflect the structure around the brightness 
peak of a non-spherically symmetric cluster. We fit a 
scaling relation of the form 

Pi 



0<xTZIS(1 + z) 



(13) 

using the orthogonal fit (Eq. |11[ ) appropriately modified 
for the additional parameter. 

Table [5] gives the Tx dependence and intrinsic scat- 
ter for seven relations using non-core subtracted quanti- 
ties (e.g., Lx~Tx), for the same relations using core sub- 
tracted quantities (e.g., LxcsrTx.es), and for the same 
relations adding the third parameter Jo (e.g., Lx-Tx-Iq)- 
The 3-paramctcr Lx relations have even lower intrinsic 
scatter than the core subtracted relations; for the M g 
relations, the reverse is true. The scatter is little differ- 
ent between the different methods for the isophotal size 
relations, with perhaps slightly lower scatter in the core 
subtracted relations. 

Interestingly, the slopes for the 3-parameter Lx rela- 
tions are even lower than those of the core subtracted 
relations, and are ~ 2er lower than the self- similar ex- 
pectation a = 2. For the M g and Ri relations, however, 
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TABLE 5 

Scatter and Slope Comparisons, Core Subtracted vs. 3-Parameter 





Standard Relation 


Core Subtracted 




3-Parainctcr 




Relation 


a 


Cint 


a 


<7int 


a 


P 




-kx,2500~Tx 


9 7 r+0.29 
„ ,r+0.33 

t.oo_ 24 
1.82 ± 0.08 


60+ 08 
39+ 12 
0.14 ± 0.02 


? 00+ 023 
z.uu_ .i 9 

2 26+ 29 

z ' ZD -0.33 

1 63+ 09 


0.28 ± 0.05 
21+ 08 

u - zl -0.07 

00+ 05 

U.UU_ 00 


1.72 ± 0.13 
1.72 ± 0.13 
1.70 ± 0.07 


0.39 ± 0.03 

n r, H + 0.06 

0.10 ± 0.02 


0.21 ± 0.03 

i>;+ 002 
0.07 ± 0.02 




1.74 ± 0.09 


0.13 ± 0.02 


1.56 ± 0.10 


0.09 ± 0.04 


1 70+ 010 


0.02 ± 0.02 


0.12 ± 0.02 




0.81 ± 0.07 
0.81 ± 0.09 


0.13 ± 0.01 
0.16 ± 0.02 


70+0- 07 

°- 66 -0.08 


0.13 ± 0.02 
0.16 ± 0.02 


0.89 ± 0.09 
0.82 ± 0.09 


0.01 ± 0.02 
-0 02+ - 03 


0.13 ± 0.02 
0.16 ± 0.02 






0.14 ±0.02 


0.65 ± 0.10 


12 +003 


0.76 ± 0.08 


-0.03 ± 0.02 


0.13 ± 0.02 



Note. — Scatter is given in base e. 

the 3-parameter slopes are consistent with those of the 
original relation, i.e., steeper than the core subtracted 
relations. 

The M g and i?j results by themselves would suggest 
that the three-parameter fit does not remove cool core- 
induced average temperature biases as completely as us- 
ing core subtracted temperatures does; i.e., the bright- 
ness of a cluster's core is not a perfect indicator of the 
strength of the core. The reduced scatter in the three- 
parameter Lx relations compared with the core sub- 
tracted relations, however, indicates that differences in 
cool core and non-cool core clusters persist outside the 
0.2r5oo core exclusion radius. Together, these results may 
lend some additional weight to the argument that cool 
core and non-cool core clusters differ in ways other than 
their apparent relaxation as determined by the develop- 
ment of a cool, dense core. 
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Fig. 15. — Central surface brightness Iq versus redshift. The 
values of Iq have been scaled by E(z)~ 3 , as described in the text, 
and are given in cgs units, i.e., erg s — 1 cm -2 arcmin - 2 . 



8.2. Evolution of Scatter 
As mentioned in § |4.2| we see a quali t ative decrease 



in scatter at higher redshifts. Kay ct al. (2007) found a 



decrease of a factor of ~3 in the luminosity-temperature 
relation in simulations, which they ascribe to merger ef- 
fects at lower redshift. However, observational studies 
have found that clusters are in fact more structurally 



disturbed at higher redshift (e.g., Jeltema et al. 20051 



Further more, c ool cores are nearly ubiquitous in the 



to observational results that find a fairly constant coo l 
core fraction of ~50 % up to z — 0.4 ( |Bauer et aTT]|2005| ); 
O'Hara et al. ( 2006 1 showed that cool core-related effects, 
and not mergers, are the primary contributors to scaling 
relation scatter at low redshift, and so clearly accurate 
simulation of core evolution is required if simulations are 
to constrain the evolution of this scatter. 

One way of gauging the effects of cool core development 
on scaling relation scatter is to look at the evolution of 
the central surface brightn ess I q. In Figure 15 we plot Iq, 
measured as described in § |8.1[ redshift. Like other clus- 
ter observables, Iq should evolve with redshift as clusters 
grow and the average density drops with the cosmic ex- 
pansion. Because Iq is a measurement of the emission 
from a cluster along the line of sight through its center, 



Iq OC 



dr, 



(14) 



and density depends on redshift as as E(z) 2 , and cluster 
radius depends on redshift as -E(z) -1 , we expect Iq cx 
E(z) 3 if clusters evolve self-similarly. Thus the values 
of Iq in Figure 
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Kay et al. (2007) simulations at all redshifts, in contrast 



are scaled by E(z)~ 3 , and if clusters 
evolve self-similarly we would expect no average change 
with redshift in I E(z)~ 3 as plotted. 

Qualitatively, however, it appears that the clusters 
with the highest Iq appear at low redshift, indicating 
a change in core structure at these redshifts. This is 
consistent with our findings that scaling relations with 
core subtracted quantities evolve faster with redshift 
than those with non-core subtracted quantities. Further- 
more, the overall scatter appears to increase at lower red- 
shifts, consistent with what we have found in observable- 
temperature relations, indicating a wider range of core 
and other structural variations as clusters develop. To- 
gether, these trends can be explained by an increasing 
cool core fraction, or an increase the the strengths of 
cool cores in those clusters that have them, at lower red- 
shifts. A larger sample of clusters would enable a more 
definitive investigation. 

9. DISCUSSION 

Our study indicates that cluster evolution is inconsis- 
tent with the simple self-similar model of cluster for- 
mation via gravitational collapse with no other heat- 
ing or cooling processes. There is a substantial body 
of observational work in this area already, so in this sec- 
tion we discuss the similarities and differences between 
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our work and earlier studies of scaling relation evolu- 
tion. The ultimate goal of such observations is to con- 
strain models of cluster formation; predictions of how 
cluster evolution will be modified by non-gravitational 
processes can be m ade both via simple analytical mod- 
els (e.g., Voit 2005 1 and from detailed h ydrodynamical 
" Muanwong et alT]|2006 1 . 



simulations (e.g 



9.1. Luminosity-Temperature 

The X-ray luminosity-temperature relation is by far 
the most studied cluster scaling relation, with several 
studies using Chandra or XMM. These studies have gen- 
erally found evolution in Lx~Tx relations that is ei- 
ther consistent with the self-similar expect a tion (e.g., 



Vikhlinin et al | [20021 |Lumb et al | pOl) |Kotov fc 



Vikhlinin 2005; Maugh an et al.l|2006D or more negative 



e.g., |Ettori it al.||20046Tpr anchcsi et al. 2007). An in 



teresting exception is |Mbrandi et al.| ( 2"007p , who found 
positive evolution when using their entire 24 cluster sam- 
ple, but marginally negative evolution when using only 
the 11 clusters which were identified as having cool cores. 

Qualitative examinatio n of the redshift scaling in our 
sample (Figures [6] [9~J and 12 1 clearly indicates the need 



to include clusters at redshifts as high as possible. Of the 
other studies mentioned above, the only ones tha t extend 



to r edshifts beyond z = .8 are Ettori et al. (2004b I 



and Branchesi et al. (20071, who find n egative evol ution 
with r espect to sell-similar, as we do; [Vikhlinin et al.| 
( 2002 1 , who see no evolution with respect to self-similar, 
but whose methods (e.g., measurement of luminosities 
within fixed 2 Mpc apertures) are quite different from 
later studies , making comparison difficult; and Maughan 



et al. (2006[), whose result is only marginally consistent 



with the self-simil ar expectation. 



The work of |Branchesi et al. ( 2007 1 in particular is 
interesting to compare to ours, because they use a Chan- 
dra sample covering a similar redshift range (though with 
only 17 members), and study two scenarios similar to 
our self-similar and no evolution scenarios. They find 
negative evolution with respect to self-similar, though 
at lower significance than our result; with an additional 
22 clusters from three other Chandra studies, the signifi- 
cance increases. In a no evolution scenario, they find the 
Lx~Tx relation evolution to be consistent with zero, as 
we do in our Lx.25oo~ Tx relation, which is mo st directly 



comparable. However, Branchesi et al. (2007) addition- 



ally measure scaling with respect to slopes and normal- 
izations from local relations, obtain poor fits, and con- 
clude that there is different evolution in the luminosity- 
temperature relation between < z < 0.3 and above this 
range. As discussed, however, there are systematic dif- 
ferences between cluster parameters measured with dif- 
ferent instruments, or even the s ame instrument in dif- 
ferent studies, as is shown in the Branchesi et al. (20071 
results where fits worsen as additional clusters are added 
from other Chandra studies. If there is a sharp change at 
low redshift, quantifying it will require a homogeneously 
reduced sample, a task made unfortunately difficult for 
Chandra by its small field of view. 

Results from simulations suggest possible explanations 
for the slower than self-similar evolution that we ob- 
serve in the Lx~Tx relation. W hile not trying to e xactly 



produced simulations using different models for the in 



duplicate observed relations, |Muanwong et al. ([2006 1 

the 



crease in entropy of the ICM. Their results show that, 
as naively expected, a simple radiative cooling model re- 
sults in faster than self-similar evolution in luminosity- 
temperature because of reduced mean cluster tempera- 
tures and increased luminosities. They found slower than 
self-similar evolution using simple preheating and stellar 
feedback models, with the latter's negative evolution sig- 
nificantly greater than the former. While their models 
are simple and cannot be directly used to test specific 
realistic models, these results do illustrate the usefulness 
o f scaling relat i ons in const raining cluster ph ysics. 



Ettori et al | ( |2004a[ ) and |Kay et al] (|2007j have stud- 
ied scaling relation evolution in simulations that include 
radiative cooling, star formation, and feedback. Both 
studies found significant (3> 3er) negative evolution with 
respect to self-simil ar in bolome tric £x,500 - ^x relations; 
specifically, |Ettori et al.] ( |2004a[ ) found 7 = -0.76 ± 0.08 
(depending on the exact method used; the other p ossi 



ble values a re the same within the uncertainty) , and Kay 
eIaT1 ( [2007| found 7 = -0.98 ±0.03 when using non-core 
subtracted quantities, and 7 = —0.61 ± 0.04 when mea- 
suring luminosities and temperatures excluding the cen- 
tral 50 kpc. Though direct comparisons may not be pos- 
sible given differences in measurement of cluster temper- 
atures between simulation and observation, differences 
in how the luminosities are measured, and the fact that 
our £x,5oo sam pl es are relatively small and consequently 
have large uncertainties in their fit parameters, the sim- 
ulation results are consistent with our results in Table [2] 
for -Lx,500 an d -kx,2500 relations. The more ne gative scal- 
i ng in the non-core subtracted relation that Kay et aL] 
(20071 found in simulations is matched by our data, and 
indicates that the primary source of the slower than self- 
similar evolution in the Lx~Tx relation is due to clusters 
being undcrluminous at higher redshifts, and not to tem- 
perature biases from cores. This slower than expected 
increase in luminosity at high redshifts indicates a po- 
tential source of difficulty for X-ray cosmology surveys, 
in that it may be more difficult to find large numbers of 
high-redshift clusters than has generally been assumed. 

9.2. ICM Mass-Temperature 

The ICM mass-temperature relation is less well stud- 
ied than luminosity-temperatu re, and results are more 
varied. | Vikhlinin et al. (2002) found significantly pos- 
itive evolution relative to the self-similar expectation 
when measuring masses within a radius defined in terms 
of the averag e baryon density of the Universe; Ettori 
et al. (2004b I found marginally significant (1-2 a) neg- 
ative evolution with respect to self-simi lar (7 = —(0.1- 
0.4), d epending on the method used); |Maughan et al.| 
( |2006[ ) claim consistency of their high-redshift sample 
with low-redshift clusters when self-similar scaling is ap- 
plied, though the y do not attempt to directly measure 
any evolution; and Morandi et al. ( |2007[ ) find significantly 
positive e volution with respec t to sell-similar. The simu- 
lations of |Ettori et al. ( 2004a I predict negative evolution 
(7 = —(0.1-0.2), depending on the method) at the 1-2 
a level. To this we compare our results, in which we 
find that M g has negative evolution with respect to self- 
similar at the 1-3 a level, depending on the radius and 
whether core subtracted parameters are used. 

9.3. Gas Fraction 
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An unchanging gas mass fraction, or one that changes 
in easily quantifiable ways, is an essential component of 
cosmological studies that use measurem ents of cluster gas 



Ettori et al. 


2003 


Allen et al. 


2004, 



varie s by cluster mass and by radius within a cluster 



(e.g. jDavid et al.|1995||Mohr et all l999 
120031 ISadat et al | |2005jr B= 

fa at hi 
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Sadat et al. 1 2TKJ5p claimed to 
ind a decrease in f g at higher redshifts when assuming a 
standard ACDM cosmology, consistent with our findings 
that f g within r25oo decreases with redshift relative to the 
self-similar expectation. The angular diameter distance 
of clusters, which is used in these cosmological studies, 

varies with f g as cLa oc fg^, and so our observed ~25% 
decrease in f g between redshifts and 1 corresponds to a 
bias of ~17% in eU estimates based on constant f g over 
the same redshift range. 

Simulations that include radiative cooling, star for- 
mation, and feedback processes likewise predict this de- 
crease in f g with redshift, with the magnitude of that 
decrease being larger at smaller fractions of the clus 



ter vi rial radius (e.g., Kravtsov et al. 2005 



Ettori et al. 



2006 1 . As with the Lx and M g evolution, the predicted 



magnitude of this evolution differs according to t he sim- 



ulation par ameters and the numerical codes used ( Ettori 
et al.|2006 l, and so observational results such as ours will 



provide constraints as simulation quality improves. 

As h as been demonstrated IFerramacho fc Blanchardl 
(20071, the results obtained from cosmological studies 
that assume constant gas fraction depend heavily on the 
radius within which measurements are made, with radii 
closer to the virial radius giving results that disagree 
greatly with the concordance model. Though measure- 
ments at large radii require extrapolation that may intro- 
duce additional biases, such results when combined with 
evidence of the radial and redshift dependence of f g give 
strong warning against ready acceptance of cosmologi- 
cal results that assume constant f g , particularly when 
measurements are made at small radii such as r 2 5oo- 

Though our results suggest difficulties for cosmologi- 
cal stud ies that assume constant f g , we have presented 
in §7.2| an alternative method for studying cosmology 
that involves using information about the evolution of 
fg to measure angular diameter distances using isophotal 
sizes. This cosmic consistency test requires joint analysis 
of cluster structure, using L x and M g to constrain ICM 
evolution and Ri to estimate distances. Consistent in- 
put and output cosmological models are guaranteed only 
around the correct model. 

10. CONCLUSIONS 

We study the evolution of the ICM using X-ray scal- 
ing relations measured from a large, homogeneously an- 
alyzed sample of clusters spanning 0.2 < z < 1.2. We 
use luminosity- and ICM mass-temperature relations, in- 
cluding both relations with and without core subtracted 
quantities, to test scenarios of standard "self-similar evo- 
lution" and of "no evolution". We also study the evo- 
lution of isophotal size-temperature relations, for which 
(under certain assumptions) these two scenarios are iden- 
tical. Finally, we compare the scatter in scaling relations 
after attempting to reduce cool core-induced scatter in 
two different ways. Our principal results appear below: 



1. Luminosity- and ICM mass-temperature relations 
evolve less rapidly than expected in the self-similar 
evolution scenario; that is, clusters at higher red- 
shifts have systematically lower luminosity and 
ICM mass at a given temperature than would be 
expected if clusters evolved self-similarly. The core 
subtracted relations have a combined consistency 
with the self-similar prediction of <0.1%; non- 
core subtracted relations are even more inconsis- 
tent with the self-similar prediction. 

2. The data are also inconsistent with the no evolu- 
tion scenario, though not at as strongly as in the 
self-similar scenario. The core subtracted relations 
evolve more rapidly than expected at higher red- 
shift in this scenario, with combined probability of 
consistency with no evolution of 1%. 

3. The evolution in the £xcs - Ixcs and M g -Txcs re- 
lations is consistent with a simple evolution in gas 
fraction, with evolution in f g at > 99% confidence 
(jf — —0.39 ± 0.13) in the self-similar evolution 
scenario when using core subtracted observables 
measured within r25oo- 

4. Isophotal size evolves with redshift at a rate that 
depends on the isophote used, reflecting evolution 
in the ICM spatial distribution in clusters. Evolu- 
tion of isophotal size at a low isophote (i.e., well 
away from the core) is consistent with that ex- 
pected given the measured f g evolution. 

5. Relations with core subtracted quantities in general 
have more positive evolution than relations with 
the cores included, suggesting that either the cool 
core fraction decreases with increasing redshift, or 
that the cool core fraction remains constant but the 
cores that do exist are weaker at high redshift. This 
is supported by direct observations of the redshift 
dependence of central surface brightness, a good 
indicator of cool core development; the scatter and 
magnitude of Iq increase at low redshift. 

6. Core subtracted relations generally have tempera- 
ture dependences that are shallower than non-core 
subtracted relations, and thus are more consistent 
with the slopes predicted by the self-similar model 
for each scaling relation. 

7. The use of core subtracted quantities for scaling 
relations and the use of non-core subtracted quan- 
tities with the addition of a third parameter, the 
central surface brightness, both significantly reduce 
scaling relation scatter by compensating to some 
extent for cool core-related effects. 

8. Scatter in observables at fixed temperature appears 
to decrease with redshift. This could indicate an 
increase in the cool core fraction, an increase in the 
strength of cool cores in those clusters that have 
them, or both. 

Cluster simulations are still improving with regard to 
their ability to accurately model non-gravitational pro- 
cesses and thus to directly test specific models by com- 
parison to observational data. However, our results of 
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negative evolution with respect to self-similar expecta- 
tions in Lx and M g , and consequently in f g , provide 
important constraints for future computational studies. 
Our findings provide new warnings with regard to the as- 
sumptions made when using f g measurements to study 
cosmology. It has long been established that f g varies 
with radius inside clusters and varies with cluster mass 
when measu red within r 500 (e.g., David et al.|1993 Mohr 



et al. 19991. Our results strongly suggest that f g varies 
with redsmft as well. Given the differences in behavior of 
collisionless dark matter and the ICM (particularly the 
ICM's sensitivity to radiative cooling and feedback from 
AGN and supernovae), perhaps it should not be surpris- 
ing that these components vary differentially with radius, 
cluster mass, and even redshift. 

At the same time, the combination of isophotal size 
measurements with measurements of the evolution of f g 
from Lx and M g relations provides a promising tool for 



measuring angular diameter distances. Our proposed 
cosmic consistency test would allow one to use cluster 
structure and its evolution to constrain cosmology in a 
manner complementary to more established techniques. 
Finally, our results underscore the need to directly cali- 
brate (or self-calibrate) mass-observable scaling relations 
in large cluster survey cosmology experiments. Cluster 
structural evolution is subject to a wide range of interest- 
ing physics, and determining that mix reliably enough for 
even the most sophisticated simulations to precisely pre- 
dict cluster mass-observable scaling relations and their 
evolution will remain enormously challenging for the fore- 
seeable future. 
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Observation and Spectral Fitting Information 
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g+i.'i 

— 5 


7-i + „i 

— 0.8 


RX J2129. 7 + 0005 . . . 





235 


552 


9.9 


21 


29 


40 


1 


+ 00 


05 


18 


218 


5.7 + 0.3 


6. 


7+1.1 
— 0.6 


<3- 8+ „'« 
— 0.6 


RX J0439. 0+0715 . . . 





245 


3583 


19.2 


04 


39 


00 


.8 


+ 07: 


15 


:58 


243 


7.4 + 0.6 


6. 


7+1.0 
— 7 


7.0+J' 3 
— 1.0 




















-10 


14 


:21 


360 




5. 


4 +0.5 


5 4+°' 4 
8 -0.3 


A1835 





252 


495 


18.4 


14 


01 


01 


9 


+ 02 


52 


41 


187 


8.2 + 0.2 


16 


.3+% 3 


16.1 + 3 " 
— 2.9 






255 






00 








+ 09: 


09 


27 


260 


8 6 + 1 - 4 
B.O_ g 


8. 


4 +l 2 9 5 
— 1.6 


8 0+ 2 ' 
8.U_ 1 g 


MS 1455.0 + 2232 .... 





258 


4192 


91.6 


14 


:57 


15 


.1 


+ 22: 


20 


34 


148 


4.7 + 0.1 


5 


6±0.3 


5.6 + 0.3 


MS 1006.0 + 1202 .... 





261 


925 


15.4 


10 


:08 


47 


.5 


+ 11: 


47 


34 


234 


6.1+0.4 


6. 


6 + 1.1 
"-0.7 


« « + l-4 
6£i -0.7 


A697 





282 


4217 


19.5 


08 


:42: 


57 


.6 


+ 36: 


21 


55 


276 


10 5 + - 9 
±uo — 0.5 


11 


.9+1.2 


11.6 + 1.3 


A611 





288 


3194 


24.3 


08 


00 


56 


.8 


+ 36: 


03 


23 


172 


8 9+ - 7 


11 


+ 3.6 
• 8 -2.2 


i2.5l 8 ;| 


ZwCl 3146 





291 


909 


43.7 


10 


23 


39 


6 


+ 04: 


11 


10 


246 


6.5 + 0.1 


8. 


7+0.7 
-0.4 


8 ^°0'.5 


A781 





298 


534 


9.9 


09 


:20 


21 


.6 


+ 30: 


30 


20 


264 


. , + 0.6 
— 0.4 


5. 


,+0.7 
d -0.4 


5 2 + 0.6 
5 -0 4 


MS 1008.1-1224 





301 


926 


28.6 


10 


10 


32 


2 


-12 


39 


23 


196 


6.4 + 0.4 


6. 


,+0.9 
°-0.6 


6 ' 6 -0.6 


RXC J2245. 0+2637 . 





304 


3287 


14.6 


22 


45 


04 


9 


+ 26: 


38 


02 


150 


5.9 + 0.3 


7. 


,+1.2 
-0.9 


6 7 +1.3 
— 0.8 


A1300 





.308 


3276 


13.7 


11 


31 


55 


.3 


-19 


54 


:46 


268 


= = + 0.7 
-0.6 


9. 


4 + 1.0 
-0.9 


9 1 + 10 
— 0.9 


A2744 





308 


2212 


22.1 


00 


1-1 


15 


3 


-30 


22 


50 


235 


10.1 + 0.6 


9. 


,+0.7 
-0.6 


9.3 + 0.7 


MS 2137.3 2353 . 





313 


5250 


25.6 


21 


40 


15 


2 


-23 


39 


38 


148 


5.0 + 0.2 


5 


+ 0.5 


5 2 + °- 3 
5 -0.6 


A1995 





318 


906 


10.0 


14 


52 


58 


.6 


+ 58: 


02 


:58 


191 


8 1+ 1 - 
S -0.8 


6. 


+ 1.0 
-0.8 


5 7+ 1 ' 1 


ZwCl 1358+6245 





327 


516 


20.0 


13 


59 


51 


.4 


+ 62: 


30 


53 


185 


9 1+0.9 
-0.8 








A1722 





328 


3278 


14.6 


13 


20 


08 


3 


+ 70 


04 


34 


203 


9 1 + 1 ' 5 
— 1 2 


13 


,+6.4 
-4.2 


10 6 + 7 ' 4 
■ —2 9 


RXC .70404. 6 + 1109 . 





355 


3269 


21.8 


04 


0-1 


33 


.7 


+ 11 


08 


25 


321 


5 6+ ' 8 
,J -0.7 


S. 


. +0.9 
-0.6 


5 1 + 10 
01 — 0.6 


RX .71532.9 + 3021 . . . 





362 


1649 


8.1 


15 


32 


54 


.0 


+ 30: 


21 


:04 


128 


6.1+0.3 


8. 


1+1.6 
1 -1.2 


7.5+i : 5 


A370 





373 


515 


53.9 


02 


39 


54 


.5 


-01 


34 


:47 


184 


8 7 + °- 5 


8 


1 + 0.5 


7.8 + 0.5 


ZwCl 1953 





374 


1959 


21.0 


08 


50 


08 


-1 


+ 36: 


04 


35 


214 


7.6 + 0.5 


6. 


+ 0.6 
°-0.5 


6.2 + 0.5 


RXC .70949.8 + 1707 . 





383 


3274 


14.3 


09 


49 


52 


.4 


+ 17: 


07 


10 


246 


7 8 + °- 7 

— 0.6 


8. 


1 +1.2 
'-1.1 


7.5 + 1.2 


C1G .71416 + 4446 .... 





400 


541 


29.9 


14 


16 


28 


4 


+ 44 


46 


42 


128 


3.8 + 0.3 


4. 


,+0.7 


4 3 + °- 8 
-0.5 


RXC .72228.6+2036 . 





412 


3285 


19.8 


22 


28 


32 


.1 


+ 20: 


37 


23 


244 


8.1 + 0.5 


7. 


„ + 0.8 
a -0.7 


8 4+ 1 ' 4 
8l -0.8 


MS 0302.7+1658 





426 


525 


10.0 


03 


05 


31 


7 


+ 17: 


10 


05 


82 


3 6 + 05 
— 4 


2. 


g+0.7 
-0.5 


2 7+°- 5 
— 4 


MS 1621.5 + 2640 .... 





426 


546 


30.0 


16 


23 


35 


.0 


+ 26 


34 


26 


197 


6 4+ - 6 
-0.5 


6. 


4 + 0.8 
-0.7 


3+0.8 
-0.7 


MACS .70417.5-1154 





440 


3270 


11.9 


04 


17 


33 


.5 


-11 


53 


58 


270 


9.4 + 0.7 


11 


4+1.9 
-1.6 


lo.ej; 2 ; 3 


RXC .71206.2-0848 . . 





.1-10 


3277 


23.4 


12 


06 


12 


.2 


-08 


48 


:05 


236 


11.4 + 0.9 


12 


k + 1 ■ 7 
' 5 -1.5 


12 7+ 2 ' 1 
— 1 8 


C1G .70329 0212 





450 


6108 


39.5 


03 


29 


41 


6 


-02 


:11 


:46 


127 


5.9 + 0.2 


6. 


o + l.l 

8 -0.6 


7 2+ 1 ' 


RX ,71347.5-1145. . . . 





451 


3592 


57.7 


13 


47 


30 


.7 


-11 


45 


11 


167 


13 4 + 0-5 
-0.3 


13 


g + 1.7 
-0.9 


12.8+1' 3 


C7G ,71701 + 6414 





453 


547 


48.2 


17 


01 


24 


.0 


+ 64: 


14 


11 


108 


4.7 + 0.3 


5. 


,+0.9 
-0.5 


5 l+°' 8 

6 -0.6 


3C 295 





461 


2254 


79.8 


14 


11 


20 


2 


+ 52: 


12 


08 


128 


5.7 + 0.2 


5. 


4 + 0.6 
-0.5 




C1G .71621 + 3810 





461 


6172 


29.8 


16 


21 


25 


.0 


+ 38: 


10 


07 


118 


™i8:2 


7. 


4 + 1.4 
-1.3 


8.2+ 2 ; 4 


C1G .71524+0957 





516 


1664 


50.1 


15 


2-1 


39 


.8 


+ 09: 


57 


46 


112 


4.8 + 0.4 


4. 


6 + 0.6 
°-0.5 


4 9+°' 7 
^■ y _0.6 


MS 0451.6-0305 





539 


902 


32.3 


04 


54: 


11 


.9 


-03 


00 


:56 


147 


9.7 + 0.8 


8. 


, + 1-1 
-0.8 


8 3+ 1 - 4 
8 ' 3 -1.0 


MS 0015.9 + 1609 .... 





541 


520 


67.4 


00 


18 


33 


.7 


+ 16 


26 


17 


197 


9.7 + 0.5 


9. 


„+0.7 
-0.6 


10 l+°' 9 
lul - 0.8 


C1G .71149+2223 .... 





544 


3589 


20.0 


11 


49: 


35 


7 


+ 22: 


24 


:04 


177 


9.8 + 0.8 


9. 


i +1.0 
1 -0.9 


9 0+ 1 ' 2 

y.u_ 9 


C1G J1423+2404. . . . 





545 


4195 


103.6 


14 


23 


47 


8 


+ 24: 


04 


41 


156 


5 4 + °- 2 
-0.1 


5 


+ 0.3 


4.6 + 0.3 


C1G J1354-0221 





546 


5835 


37.5 


13 


54 


17 


2 


-02 


21 


50 


94 


4 j+0.8 
* -0.3 


4. 


„+l.l 
-0.9 


3 9+ 12 
3J -1.0 


C1G .70717+3745 .... 





548 


4200 


59.1 


07 


17 


31 


.3 


+ 37: 


45 


35 


244 


11 5+°' 7 
11 -0.8 


10 


, + 0.8 
" 5 -0.6 


9.9 + 0.6 


C1G .71120+2326 .... 





562 


1660 


69.3 


11 


20 


57 


.5 


+ 23: 


26 


:34 


128 


4.2+0-6 


4 


7 + 0.7 


3.9 + 0.4 


C1G .72129-0741 





570 


3595 


19.9 


21 


29 


26 


2 


-07 


41 


28 


166 


11.8+" 


9. 


0+2-7 
"-1.2 




MS 2053.7 0449 





583 


1667 


44.5 


20 


56 


21 


3 


-04 


37 


:49 


69 


4 0+ - 5 
-0.2 


3. 


6 + 0.8 
-0.5 


3.2+0; 8 


C1G J0647+7015 .... 





584 


3584 


19.9 


06 


:47: 


50 


.6 


+ 70: 


14 


:54 


160 


i..o+3;» 








C1G .70542-4100 





634 


914 


48.6 


05 


42 


49 


.6 


-40 


59 


58 


118 




5. 


4 + 1.0 
-0.6 


6 2+ 12 
°-- ! — 1.0 


C7G .71419+5326 





640 


3240 


9.1 


14 


19 


12 


2 


+ 53: 


26 


09 


59 


4 l+°- 8 
-0.6 


3. 


4 + 0-8 
-0.7 


3 1 + 1 - 6 


C1G .70744+3927 





686 


6111 


49.5 


07 


1-1 


52 


.8 


+ 39: 


27 


27 


118 


9.6 + 0.9 


11 


7+2.2 
-2.0 


10 4+ 3 ' 2 


C7G .71221 + 4918 .... 





700 


1662 


78.3 


12 


21 


25 


9 


+ 49: 


18 


28 


138 


6 5 + °- 8 
°' 5 -0.6 


6. 


4 + 1-0 
-0.7 


6 1+ ' 9 


C1G .71113 2615 





730 


915 


62.5 


11 


13 


05 


.0 


-26 


15 


:40 


79 


, 7 +0.6 


2. 


„+0.6 
8 -0.4 


2 -6i°0'.4 


C1G 1137+6625 . . 





.782 


536 


27.6 


11 


:40 


22 


.4 


+ 66: 


08 


16 


98 


5 9+ 1 ' 2 
°' 9 -0.9 


6. 


,+2.5 
1 -1.7 


6. 4 i2'.2 


RX .71350.0 + 6007. . . 





804 


2229 


58.3 


13 


:50 


48 


.3 


+ 60: 


07 


11 


98 


4 l+°- 8 
-0.6 


4. 


„+1.6 
-0.8 


4 5+ 2 'l 
— 1.2 


RX ,71317+2911 .... 





805 


2228 


111.3 


13 


17 


21 


.8 


+ 29: 


11 


19 


69 


3 s+1-7 
3 ' 8 -0.9 


3. 


, + 3.1 


, , + 3.0 
2 - -0.5 


RX .71716 + 6708 





813 


548 


51.2 


17 


16 


49 


.1 


+ 67: 


08 


2-1 


108 


6 4+°' 9 
°' 4 -0.8 




6 + 1.2 
"-0.8 


6 4+ 2 ' 3 


C1G .71056-0337 





.826 


512 


66.7 


10 


56 


59 


.5 


-03 


37 


34 


118 


9 2 + 1.5 
1 -1.2 


8. 


7+1.8 
-1.1 


a 6+ 2 - 3 
8 ' 6 -1.5 


C1G .71226+3332 .... 





890 


3180 


31.6 


12 


26 


58 


.0 


+ 33: 


32 


46 


108 


12.2 + J- 8 


13 


+ 4.0 
'°-3.2 


10.3+3' 1 


C1G .71415 + 3611 


1 


030 


4163 


89.2 


14 


15 


11 


.2 


+ 36: 


12 


03 


79 


6 S+ 1 ' 
6 ' 8 -0.7 


6. 


, + 1.8 
2 -l.l 


— 1.6 


C1G .71252 2927 


1 


235 


4198 


162.5 


12 


52 


54 


.4 


-29 


27 


16 


69 


5 7+1- 4 
5 -1.0 


5. 


, + 1.6 
• -1.0 


5 2 + 2 ' 4 
6 ' -1.3 



a Exposure time after light curve filtering.^ 3 Coordinates given arc center of spectral extraction aperture. c Core-subtracted temperature 
measured assuming self-similar evolution of rg. d Core-subtracted temperature measured assuming no evolution of r$ . 
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TABLE 7 
(3 Model Parameters 



Cluster Fit aperture 1^ R c i 1^ R c 2 

(arcscc) (org s " 1 cm " 2 arcmin " 2 ) (arcs'cc) (erg a " 1 cm " 2 arcmin " 2 ) (arcsUc) 



A665 


394 


0.62 + 0.01 


1.9 + 0.0 E-12 


65.6+i;* 






A963 


197 


0.55 + 0.00 


6.6 + 0.1 E-12 


2i.i+0;| 






RX .10439.0+0520. . . 


148 


67+ - 04 
UD -0.02 


3.0+°|E-12 


28.0+11 


5.4 + 0.4 E-ll 


5 4+ ' 5 
' 5 - -0.4 


A1423 


187 


0.46 + 0.01 


7.2 + n'cE-12 
— 0.5 


10.5±0.8 






ZwCl 2701 


153 


0.58 + 0.01 


1.5 + 0.0 E-ll 


12.3±0.3 






A773 


256 


0.60 + 0.01 


2.6 + 0.1 E-12 


41.2 + 1* 






A2261 


148 


55 + - 01 
-0.00 


1.2 + 0.0 E-ll 


18.1+°? 






ACO 2246 




0.52 + 0.01 


3.7 + 0.3 E- 1 2 


q n + O.S 






A1682 


177 


0-56 + 0-06 
— 0.04 


1.0 + 0.1 E-12 


49.2 + --I 
— 7.5 






A2111 


295 


0.58 + 0.02 


1.2 + 0.1 E-12 


48.7+1;! 










62 + 0.01 


3 3 + 01 E-12 


33.3+1 2 






RX J2129. 7+0005 . . . 


157 


0.60 + 0.01 


7.4 + 0.9 E-12 


23. 4+ 2,3 
— 2.0 


6.2 + n'oE-ll 
— 0.3 


4.1 + 0.4 


RX J0439 0+0715 




61 + 0.01 


6.0 + q' 2 E-12 


26 1 + °' S 
— 1.1 






A521 


295 


0.75 + 0.00 


5.5 + 0.2 E-13 


122.0 + 2.4 






A1835 


167 


0.73 + 0.01 


5.7 + 0.2 E-12 


44.8+1;= 


1.1+0.0 E-10 


8.9 + 0.2 


A68 


246 


75 + °- D3 
U - ' -0.02 


2.3 + 0.1 E-12 


53.o±l;g 






MS 1455.0+2232 


148 


0.61+0.00 


6.4 + 0.1 E-ll 


8.9 + 0.1 






MS 1006.0 + 1202 .... 


216 


0.70 + 0.02 


1.7 + 0.1 E-12 


, n „ + 2.8 
48.0T 2 . 5 






A697 


256 


0.64 + 0.01 


3.7 + 0.1 E-12 


-1-1 6 










0.60 + 0.01 


8.4 + 0.3 E-12 


18.4 + 0.6 






ZwCl 3146 


246 


6S + - 01 
u ° -0.00 

1 47+0-36 

-0.24 


1.6 + 0.1 E-ll 


23.8 + 0.6 


1.2 + 0.0 E-10 


5.4 + 0.1 


A781 


226 


6.9 + 0.3 E-13 


157 2 + 28.0 
-21.2 






MS 1008.1-1224 


197 


0.65 + 0.02 


2.4 + 0.1 E-12 


35. Oil. 8 






RXC J2245. + 2637 . 


148 


0.66 + 0.02 


8.5 + 0.8 E-12 


21 3 + 22 
— 1.6 


4.1 + g '5E-II 


4 1+ - 6 
*■ -0.5 


A1300 


207 


0.49 + 0.01 


4-»±S:i E - 12 


22.1+1;! 






A2744 


406 


1.10 + 0.04 


1.8 + 0.0 E-12 








MS 2137.3 2353 


138 


0.64+0.01 


3.9±J;jB-ll 


10.2 + 1.1 


1.1 + 0.0 E-10 


3 3+ 05 
J ' -0.6 


A1995 


216 


0.82 + 0.03 


3.5 + 0.1 E-12 


43.5+J; 2 ! 






ZwCl 1358+6245 


157 


n ,.,.+0.03 

°- 66 -o.oi 


2.8 + 0.2 E-12 


31.4+ 24 


4 °to.4 E - 11 


3 8+ - 4 


A1722 


148 


64+ - 05 


4-0 4 
21 Io.6 E " 12 


30 7+ 7 ' 9 
iu '-3 .9 


2.9il; 2 E-12 


7 2+ 4 - 4 
'-2.5 


RXC J0404. 6+1109 . 


128 


46 + ' 04 


1.0 + 0.1 E-12 


28.6 + = -9 






RX J1532. 9+3021 . . . 


118 


0.61 + 0.01 


1.1 + 0.0 E-10 


7 8 + °' 2 






A370 


187 


0.81 + 0.02 


1.7 + 0.0 E-12 


59.4+J; 






ZwCl 1953 


246 


0.65 + 0.01 


4.5 + 0.2 E-12 


30.9+1! 






RXC J0949. 8 + 1707 . 


153 


0.63 + 0.02 


6.0 + 0.3 E-12 


27.8+lj 






C1G J1416+4446 


148 


— 0.02 


3.0 + n'^E-12 
— 0.7 


— 2.4 


2.3 + n '?E-ll 
— 0.6 


1.9+°% 
— 0.5 


RXC J2228. 6 + 2036 . 


172 


0.64 + 0.02 


5.3 + 0.2 E-12 


31.0+1' 8 , 
— 1.7 






MS 0302.7+1658 .... 


98 


0.54 + 0.02 


1 5+°-?E-ll 
— 0.2 


6.9 + 1? 
— 1.0 






MS 1621.5 + 2640.... 




67+0.03 


1.3 + 0.1 E-12 


41 6 + 3 ' 6 






MACS J0417. 5-1154 


192 


0.65 + 0.02 


4.4 + 0.3 E-12 


48.4+ 2 .f 


1.1+0.1 E-10 


5.2 + 0.3 


IV-A. 1 ^/ JlzUu.Z-Uo4o . , 




U . oUzt u .U 1 


i t -I- n i ii 


ls)1 _0.6 






C1G J0329 0212 


118 


0.52 + 0.00 


1.2 + 0.1 E-10 


3.8 + n .1 
— 0.2 






RX J 1347. 5-1145. . . . 


189 


0.65 + 0.00 


4.5 + 0.2 E-ll 


16.7+Si 
-0.5 


3.9 + 0.1 E-10 


3.7 + 0.1 


C1G J1701 + 6414 


148 


0.58 + 0.02 


i.i+.°; 2 e-12 


29.4+|l 


2.9+_°|e-11 


2.7 + 0.3 


3C 295 


118 


0.63 + 0.01 


9 ' 3 ~*"l 4 E " 12 


13 1 + 1 - 4 
-1.3 


1.3 + 0.1 E-10 


2 7+ - 2 
-0.1 


C1G J 162 1 + 38 10 . 


108 


0.60 + 0.02 


7.6 4 " 2 'gE-12 


14 S + 2 ' 5 
-2.0 


9.2+l -2 E-ll 


2 5+°- 4 
z - -0.3 


C1G J1524+0957. . . . 


128 


95 + - 14 
u y -0.10 


1.1+0.1 E-12 


56 2 + 8 - 2 






MS 0451 6 0305 




85 + 0.02 


9 2 + 2 E- 1 2 








MS 0015.9 + 1609 .... 


216 


0.70 + 0.01 


6.5 + 0.2 E-12 


37.5 + ^1 
— 1.0 






C1G J 1149+2223 .... 


295 


0.65 + 0.02 


4.3 + 0.2 E-12 


40.9+ 2 -5 
— 2.2 






C1G J1423+2404. 


98 


0.65 + 0.01 


4.2 + 0.4 E-12 


22 0+ 1 4 
— 1.0 


2.3 + 0.0 E-10 


3.7 + 0.1 


C1G J1354-0221 


157 


76 + - 12 
u -' -0.08 


O O + l'0p in 

-0.8 


39 g+S.3 
— 6.6 






C1G J0717+3745 


187 


0.82+0.02 


4.7 + 0.1 E-12 


65.6 + 2 ^ 






C1G J 1 120+2326 .... 


148 


1 74+?^ 4 
' — 0.31 


9 9+2"^ E-13 
•j ■<> 0.4 


88 

— 12.6 






C1G J2129 0741 


166 


D fi9 + 0. 02 


1.1 + 0.1 E-ll 


18.6±1.2 






MS 2053.7-0449 . 


89 


63 + ' 05 


3.9+_°;|E-12 


15.6±2i 






C1G J0647 + 7015 


148 


0.63 + 0.02 


1.3 + 0.1 E-ll 


18.4+1;3 






C1G J0542 4100 


112 


0.58 + 0.03 


2.9 + 0.2 E-12 


22 - 5 -2.4 






C1G J1419 + 5326 


118 


0.60 + 0.03 


1.9+q- 4 E-11 


7 3 + 1 ' 2 






C1G J0744+3927. . . . 


148 


0.56 + 0.01 


4.4 + 0.2 E-ll 


8.5 + 0.4 






C1G J1221 + 4918 


157 


73 + - 04 
u - '' 3 — 0.03 


2.5 + 0.1 E-12 


35.7i 2 ;l 






C1G J1113 2615 


98 


Q 73 + 0.08 
U - 7d -0.06 


4.3 + 0.5 E-12 


15.8+J-j 






C1G 1137+6625 . . . 


98 


n „ c +0.04 
°- 65 -0.03 


i.5±g;?B-ii 








RX J1350. + 6007. . . 


148 


61 + ' 05 
u - Di -0.04 


2.3 + 0.3 E-12 


21 4+ 39 






RX J1317+2911 


89 


U- ° -0.20 


5.3+^ '^E-13 


29.4+^S 


4.8+ 2 " 2 E-12 


4 3 + 2 - 4 
* d -1.6 


RX J1716 + 6708 


148 


n fiS + 0.04 
°- 68 -0.03 


8.3 + 0.6 E-12 


17.311;? 






C1G J1056-0337 . . 


197 


0.67+0.00 


5.8 + 0.3 E-12 


31.9±0.9 






C1G J1226 + 3332 .... 


128 


0.68 + 0.02 


3.3 + 0.2 E-ll 


14.5+1.0 






C1G J1415 + 3611 


98 


n 7 =+0.06 
°- 7S -0.04 


9.5 + 1.2 E-12 


18.1 + 2.4 


6 - 5+ . i' 4 B -ll 


2 5+ - 6 
25 -0.5 


C1G J1252 2927 


89 


0.54 + 0.03 


1.1+0.2 E-ll 


8.8_ 1 4 







O'Hara et al. 



TABLE 8 

Cluster Measurements Assuming Self-Similar Evolution 



Cluster 


-^X,2500 


Lx,500 




M gi50 o 


fll.BE- 


13 


R&E-14 


H3E-14 




(1O 44 L ) 


(1O 44 L ) 


(1O 13 M ) 


(1O 13 M ) 




(Mpc 


) 


(Mpc) 


(Mpc) 


A665 


3.27+0.03 




3.93+0.09 


13.05+0.30 


n 


40i0 


.02 








A963 


2.95+0.03 




2.86+0.08 


9.27+0.27 


n 


33i0 


.02 








RX J0439. 0+0520 . . 


2.17+0.02 




1.42 + 0.07 


3.79+0.19 


o 


23i0 


.01 


0.31 + 


.02 


0.40+0.04 


A1423 


1.98+0.04 




2.08+0.07 


8.01+0.28 


o 


28i0 


.01 


0.46+0 


.04 




ZwCl 2701 


2.10+0.01 




1.58+0.07 


4.63+0.21 


o 


24i0 


.01 


0.34+0 


.02 




A773 


2.92+0.03 




3.55+0.09 


11.34+0.30 





36±0 


.02 








A2261 


4.69+0.03 




3.81+0.11 


12.20+0.36 





3S±0 


.02 








ACO 2246 


0.46+0.01 




0.56+0.04 


1.91+0.15 








0.21 + 


,01 




A1682 


1.61+0.09 




2.01 + 0.08 


7.55+0.31 


o 


32±0 


.02 


0.52+0 


.05 




A2111 


1.77+0.05 




2.59+0.08 


9.03+0.29 


o 


31 iO 


01 


0.49+0 


,03 


0.66+0.06 


A267 


2.37+0.03 




2.73+0.09 


8.15+0.27 


o 


30i0 


.01 


0.43+0 


.03 


0.57+0.05 


RX .72129. 7+0005.. 


4.26+0.03 




2.75+0.11 


8.13+0.33 


o 


34i0 


.02 


0.49+0 


,04 


0.63+0.07 


RX J0439. 0+0715. . 


3.28+0.04 




3.20+0.11 


9.43+0.31 


o 


35 iO 


.02 


0.49+0 


.03 


0.61+0.05 


A521 


2.00+0.07 


3.43+0.07 


2.38+0.10 


10.80+0.44 


o 


43i0 


.03 


0.68+0 


.08 


0.82+0.11 


A1835 


10.50+0.05 




5.35+0.16 


13.56+0.42 


o 


■44±0 


02 








A68 


2.84+0.07 




3.68+0.11 


9.96+0.29 


o 


35±0 


.01 


0.48+0 


.03 




MS 1455.0+2232 . . . 


4.95+0.01 


5.45+0.04 


2.20+0.12 


6.05+0.33 


o 


30i0 


.02 


0.41+0 


.04 




MS 1006.0+1202 . . . 


1.85+0.03 




2.37+0.10 


7.21 + 0.31 


o 


29±0 


.01 


0.41 + 


.03 




A697 


5.32+0.09 




5.96+0.16 


18.56+0.50 


o 


51i0 


.02 


0.71 + 


.05 


0.88+0.07 


A611 


2.91+0.03 




3.29+0.11 


9.50+0.31 


o 


32±0 


.01 


0.44+0 


.02 




ZwCl 3146 


9.16+0.02 




4.19+0.19 


10.77+0.48 


o 


41i0 


.03 


0.56+0 


.05 




A781 


1.43+0.09 




2.05+0.12 


7.99+0.45 


o 


37i0 


.03 


0.59+0 


.07 


0.74+0.10 


MS 1008.1-1224 


1.93+0.03 




2.44+0.11 


7.44+0.35 


o 


33i0 


.02 


0.45+0 


.03 


0.56+0.05 


RXC ,72245.0+2637 


3.47+0.02 


3.85+0.05 


2.54+0.13 


6.82+0.35 


n 


30i0 


.02 


0.41+0 


.03 


0.51 + 0.04 


A1300 


4.04+0.09 


5.55+0.11 


4.28+0.15 


15.90+0.56 


o 


46i0 


02 


0.69+0 


.05 




A2744 


4.74+0.10 




6.10+0.19 


17.12+0.52 


o 


5SiO 


.03 








MS 2137.3-2353 .... 


5.15+0.02 




2.27+0.14 


5.87+0.37 


o 


29±0 


.02 


0.39+0 


,03 


0.49+0.05 


A1995 


2.88+0.04 




3.45 + 0.14 


8.43+0.33 


o 


34i0 


.01 


0.44+0 


,02 


0.56+0.04 


ZwCl 1358+6245 . . . 


2.64+0.04 




3.28+0.12 


9.21+0.33 


o 


31i0 


01 


0.43+0 


,02 


0.54+0.03 


A1722 


1.80+0.05 




2.90+0.10 


8.53+0.31 


o 


2SiO 


.01 


0.40+0 


,02 


0.54+0.03 


RXC J0404.6+1109 


1.20+0.08 


2.30+0.14 


1.73+0.11 


7.61 + 0.48 





2S±0 


.01 


0.64+0 


.07 




RX .71532.9+3021 . . 


8.84+0.03 




3.62+0.22 


9.81+0.59 


n 


36i0 


.02 


0.50+0 


.04 


0.64+0.07 


A370 


2.59+0.05 




3.84+0.16 


11.46+0.49 


n 


40i0 


.02 








ZwCl 1953 


3.43+0.06 


4.31+0.09 


3.68+0.18 


11.27+0.56 


o 


39 zbO 


.02 


0.56+0 


.04 


0.71 + 0.07 


RXC ,70949.8+1707 


4.08+0.07 


5.11+0.09 


4.05+0.20 


12.40+0.61 


o 


41i0 


.02 


0.57+0 


.04 


0.76+0.07 


C1G J1416+4446 . . . 


1.10+0.02 


1.44+0.05 


1.14+0.12 


3.78+0.40 


o 


24i0 


.01 


0.34+0 


.03 


0.44+0.05 


RXC ,72228.6+2036 


4.15+0.08 


5.66+0.08 


4.43+0.23 


13.92+0.71 


o 


50i0 


.03 


0.69+0 


.06 


0.82+0.08 


MS 0302.7+1658 . . . 


1.41+0.02 




1.15+0.14 


3.76+0.45 


o 


,23±0 


01 


0.33+0 


,03 


0.43+0.05 


MS 1621.5+2640 . . . 


1.41 + 0.05 


2.09+0.09 


2.15 + 0.14 


7.60+0.51 


o 


.32±0 


.02 


0.48+0 


,04 


0.70+0.08 


MACS J0417. 5-1154 


10.27+0.19 


13.17+0.11 


7.14+0.33 


24.15+1.13 





.63±0 


.04 


0.82+0 


.07 


0.97+0.10 


RXC J 1206. 2-0848 . 


8.02+0.10 


9.47+0.11 


7.24+0.28 


21.76+0.84 


o 


54i0 


.03 


0.73+0 


.05 


0.92+0.07 


C1G J0329-0212 


5.23+0.03 


5.92+0.09 


2.89+0.22 


9.37+0.72 


o 


37i0 


.02 


0.52+0 


,05 


0.65 + 0.07 


RX J1347. 5-1145 . . . 


18.85+0.05 




10.10+0.34 


25.94+0.87 


o 


56i0 


.02 








C1G .71701+6414... 


1.22+0.03 




1.41 + 0.14 


5.16+0.50 


n 


26i0 


.01 


0.38+0 


,03 


0.50+0.05 


3C 295 


2.91+0.01 




2.06+0.17 


5.62+0.45 


o 


27i0 


.01 


0.36+0 


,02 


0.46+0.04 


C1G J1621+3810 . . . 


3.09+0.03 




2.68+0.18 


7.99+0.54 


o 


.31±0 


.01 


0.46+0 


.03 


0.59+0.05 


C1G J1524+0957 . . . 


0.92+0.04 


1.47+0.08 


1.47+0.16 


5.28+0.57 


o 


29±0 


.02 


0.43+0 


.04 


0.52+0.06 


MS 0451.6-0305 .... 


6.09+0.09 


7.06+0.16 


6.18+0.34 


15.86+0.88 





■4S±0 


.02 


0.62+0 


.04 




MS 0015.9+1609 . . . 


5.76+0.08 


7.61+0.11 


6.14+0.34 


19.43+1.08 


o 


55i0 


.03 








C1G J1149+2223 . . . 


4.78+0.14 


7.07+0.14 


5.66+0.32 


20.17+1.12 


o 


56i0 


.03 


0.86+0 


.08 


1.08+0.12 


C1G J1423+2404 . . . 


5.54+0.02 


6.23 + 0.10 


2.75+0.28 


7.91+0.80 


o 


34i0 


.02 


0.49+0 


.05 




C1G J1354-0221 


0.53+0.03 




0.94+0.13 


3.58+0.48 


o 


20i0 


.01 


0.33+0 


.03 


0.44+0.05 


C1G J0717+3745 . . . 


7.78+0.23 


11.13+0.13 


8.51+0.41 


29.57+1.41 


o 


6SiO 


.04 


0.88+0 


.07 




C1G J1120+2326 . . . 


0.78+0.05 




1.22 + 0.16 


4.41+0.58 


o 


26i0 


.02 


0.37+0 


.03 


0.43+0.04 


C1G J2129-0741 


4.22+0.17 




5.56+0.27 


16.45+0.80 


0. 


42±0 


01 


0.63+0 


.03 


0.80+0.05 


MS 2053.7-0449 .... 


0.89+0.02 




1.13+0.17 


3.58+0.52 


o 


23i0 


.01 


0.33+0 


.03 


0.40+0.04 


C1G J0647+7015 . . . 


5.02+0.13 




7.15 + 0.28 


20.11+0.78 


o 


43i0 


.01 


0.59+0 


.02 




C1G J0542-4100 .... 


1.33+0.05 


1.95+0.08 


2.08+0.21 


7.58+0.75 





.34±0 


.02 


0.47+0 


,03 


0.58+0.05 


C1G J1419+5326 . . . 


1.40+0.04 




1.28+0.20 


3.78+0.59 





23±0 


.01 


0.33+0 


.03 


0.43+0.05 


C1G J0744+3927 . . . 


5.49+0.09 




5.03+0.36 


15.80+1.13 





,47±0 


.02 


0.64+0 


,04 


0.79+0.06 


C1G J1221+4918 . . . 


1.47+0.06 


2.18 + 0.07 


2.27+0.24 


8.27+0.89 





37±0 


.02 


0.55+0 


.05 




C1G J1113-2615 


0.64+0.02 




0.91+0.18 


2.64+0.52 


0, 


20±0 


.01 


0.27+0 


.02 


0.35+0.03 


C1G 1137+6625 


1.88 + 0.07 




2.12+0.28 


6.17+0.82 


0. 


.29±0. 


.01 


0.41+0 


.03 


0.53+0.05 


RX J1350. 0+6007. . 


0.60+0.05 


1.04+0.05 


1.01 + 0.20 


4.05+0.80 


0. 


,26±0 


.02 


0.41+0 


.04 


0.53+0.07 


RX J1317+2911 


0.20+0.02 




0.52+0.11 


1.77+0.38 








0.21 + 


.01 


0.32+0.03 


RX J1716+6708. . . . 


1.64+0.05 




2.25+0.29 


6.77+0.86 





31±0 


.02 


0.43+0 


.03 


0.54+0.05 


C1G J1056-0337. . . . 


3.23+0.23 




4.53+0.41 


16.57+1.49 


0, 


51±0 


.03 


0.62+0 


.04 


0.72+0.06 


C1G J1226+3332 . . . 


4.80+0.11 




5.94+0.44 


15.86+1.16 





,44±0. 


.02 


0.59+0 


.03 


0.71+0.04 


C1G J1415+3611 . . . 


1.48+0.03 




2.20+0.33 


6.33+0.96 


0. 


29±0 


.01 


0.40+0 


,02 


0.50+0.04 


C1G J1252-2927. . . . 


0.63+0.04 




1.28+0.28 


4.63+1.01 





.2S±0 


.01 


0.42+0 


.03 


0.54+0.05 



Evolution of the ICM 



TABLE 9 

Cluster Measurements Assuming No Evolution 



Cluster 


£x,2500 


-^g,2500 






(1O 44 L ) 


(1O 13 M ) 


(1O 13 M ) 


A665 


3.78+0.04 


4.13+0. 


09 


13.42+0.31 


A963 




2.99+0. 


09 


9.63+0.28 


RX J0439. 0+0520 . . 


2.45 + 0.03 


1.46+0. 


07 


3.84+0.19 


A1423 


2.33+0.05 


2.22+0. 


08 


8.53+0.30 


ZwCl 2701 


2.38+0.02 


1.64+0. 


07 


4.76+0.21 


A773 


3.40+0.04 


3.74+0. 


10 


11.72+0.31 


A2261 


5.41+0.03 


3.99+0. 


12 


12.70+0.37 


ACO 2246 




0.59+0. 


05 


2.01+0.16 


A1682 


1.95+0.08 


2.17+0. 


09 


7.94+0.33 


A2111 


2.11+0.06 


2.78+0. 


09 


9.43+0.30 


A267 


2.75 + 0.03 


2.86+0. 


09 


8.38+0.27 


RX J2129. 7+0005.. 


4.91+0.03 


2.87+0. 


12 


8.37+0.34 


RX J0439.0+0715 . . 


3.83+0.04 


3.35+0. 


11 


9.71+0.32 


A521 


2.58+0.09 


2.71+0. 


11 


11.42+0.47 


A1835 


12.13+0.06 


5.50+0. 


17 


13.62+0.42 


A68 


3.33+0.06 


3.85+0. 


11 


10.04+0.30 


MS 1455.0+2232 . . . 


5.74+0.01 


2.27+0. 


13 


6.20+0.34 


MS 1006.0+1202 . . . 


2.23+0.04 


2.52+0. 


11 


7.38+0.31 


A697 


6.48+0.08 


6.38+0. 


17 


19.24+0.52 


A611 


3.46+0.05 


3.45+0. 


11 


9.84+0.32 


ZwCl 3146 


10.85+0.02 


4.32+0. 


19 


10.88+0.49 


A781 


1.97+0.11 


2.39+0. 


13 


8.03+0.45 


MS 1008.1-1224 


2.41+0.04 


2.61+0. 


12 


7.69+0.36 


RXC J2245. 0+2637 


4.16+0.03 


2.65+0. 


14 


6.94+0.36 


A 1 300 


5.09+0.09 


4.71+0. 


16 


17.25+0.60 


A2744 




6.67+0. 


20 


16.77+0.51 


MS 2137.3-2353 .... 


6.15+0.03 


2.33+0. 


15 


5.98+0.37 


A1995 


3.51+0.05 


3.59+0. 


14 


8.33+0.33 


ZwCl 1358+6245 . . . 


3.25+0.05 


3.47+0. 


12 


9.46+0.34 


A1722 


2.22+0.05 


3.10+0. 


11 


8.84+0.32 


RXC J0404. 6+1109 


1.69+0.11 


2.01+0. 


13 


8.59+0.54 


RX J1532. 9+3021 . . 


10.90+0.04 


3.78+0. 


23 


10.15+0.61 


A370 


3.49+0.07 


4.25+0. 


18 


11.62+0.50 


ZwCl 1953 


4.48+0.06 


4.02+0. 


20 


11.78+0.58 


RXC J0949. 8+1707 


5.30+0.08 


4.42+0. 


22 


13.00+0.63 


C1G .71416+4446... 


1.48+0.03 


1.27+0. 


13 


4.06+0.43 


RXC J2228. 6+2036 


5.71+0.11 


4.92+0. 


25 


14.68+0.75 


MS 0302.7+1658 . . . 




1.27+0 


15 


4.09+0.49 


MS 1621.5+2640 . . . 


1.98+0.06 


2.50+0. 


17 


8.13+0.54 


MACS J0417. 5-1154 


14.27+0.22 


8.17+0. 


38 


25.86+1.21 


RXC J1206. 2-0848 . 


10.79+0.12 


7.91+0 


30 


23.15+0.89 


C1G J0329-0212 


6.94+0.04 


3.19+0. 


24 


10.32+0.79 


RX J1347. 5-1145... 


24.62+0.06 


10.58+0.35 


26.60+0.89 


C1G J1701+6414 . . . 


1.71+0.03 


1.64+0.16 


5.69+0.55 


3C 295 


3.86+0.02 


2.19+0.18 


5.85+0.47 


C1G J1621+3810... 


4.17+0.04 


2.92+0.20 


8.52+0.58 


C1G J1524+0957... 


1.48+0.05 


1.83+0.20 


5.43+0.59 


MS 0451.6-0305 .... 


8.82+0.11 


6.74+0.37 


15.49+0.86 


MS 0015.9+1609 . . . 


8.78+0.10 


7.14+0.40 


20.52+1.15 


C1G .71149+2223... 


7.55+0.18 


6.86+0.38 


22.15+1.23 


C1G J1423+2404... 


7.87+0.05 


3.03+0.31 


8.36+0.84 


C1G J1354-0221 .... 


0.85+0.05 


1.19+0.16 


3.87+0.52 


C1G J0717+3745 . . . 


12.73+0.23 


10.46+0.50 


31.03+1.48 


C1G .71120+2326 . . . 


1.35+0.05 


1.60+0. 


21 


4.17+0.55 


C1G J2129-0741 


6.34+0.21 


6.23+0. 


30 


17.66+0.86 


MS 2053.7-0449 .... 


1.37+0.03 


1.31+0 


19 


3.88+0.57 


C1G .70647+7015 . . . 


7.34+0.16 


7.89+0. 


31 


21.30+0.83 


C1G .70542-4100 .... 


2.29+0.08 


2.59+0. 


26 


8.68+0.86 


C1G J1419+5326... 


2.13+0.05 


1.44+0. 


23 


4.12+0.65 


C1G .70744+3927 . . . 


8.90+0.13 


5.86+0. 


42 


18.02+1.29 


C1G .71221+4918... 


2.69+0.07 


2.97+0. 


32 


9.09+0.97 


C1G J1113-2615 .... 


1.02+0.04 


1.06+0. 


21 


2.73+0.54 


C1G 1137+6625 .... 


3.23+0.10 


2.46+0. 


33 


6.65+0.89 


RX J1350. 0+6007. . 


1.31+0.07 


1.41+0. 


28 


4.90+0.97 


RX J1317+2911 


0.38+0.04 


0.68+0. 


15 


1.87+0.40 


RX J1716+6708. . . . 


2.96+0.08 


2.70+0. 


34 


7.26+0.92 


C1G .71056-0337.... 


6.67+0.16 


6.14+0. 


55 


19.02+1.71 


C1G J1226+3332 . . . 


8.67+0.17 


6.71+0. 


49 


16.52+1.21 


C1G .71415+3611 . . . 


3.03+0.08 


2.68+0 


41 


6.53+0.99 


C1G J1252-2927.... 




1.90+0 


41 


6.25+1.36 



